Thermal Photons and Lepton Pairs 
from Quark Gluon Plasma and Hot Hadronic Matter 

Jan-e Alam^'Q, Sourav Sarkar**, Pradip Roy'^, T. Hatsuda"'0 and Bikash Sinha' 



b,c 



a) Physics Department, Kyoto University, Kitashirakawa, Kyoto 606-8502, Japan 

b) Variable Energy Cyclotron Centre, 1/AF Bidhan Nagar, Calcutta 700 064 India 

c) Saha Institute of Nuclear Physics, 1/AF Bidhan Nagar, Calcutta 700 064 India 

Abstract 

The formulation of the real and virtual photon production rate from strongly 
interacting matter is presented in the framework of finite temperature field theory. 
The changes in the hadronic spectral function induced by temperature are discussed 
within the ambit of the Walecka type model, gauged linear and non-linear sigma 
models, hidden local symmetry approach and QCD sum rule approach. Possibility 
of observing the direct thermal photon and lepton pair from quark gluon plasma 
has been contrasted with those from hot hadronic matter with and without medium 
effects for various mass variation scenarios. At SPS energies, in-medium effects 
of different magnitude on the hadronic properties for the Walecka model, Brown- 
Rho scaling and Nambu scaling scenarios are conspicuously visible through the low 
invariant mass distribution of dilepton and transverse momentum spectra of photon. 
However, at RHIC energies the thermal photon (dilepton) spectra originating from 
Quark Gluon Plasma overshines those from hadronic matter for large transverse 
momentum (invariant mass) irrespective of the models used for evaluating the finite 
temperature effects on the hadronic properties. It is thus expected that both at 
RHIC and LHC energies the formation of Quark Gluon Plasma in the initial stages 
may indeed turn out to be a reahstic scenario. 
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In the following we list some of the symbols which have been used. 

Afj^i, transverse projection tensor 
B^i, longitudinal projection tensor 

free vacuum propagator for spin 1 particles 
D^^ matrix propagator for spin 1 particles at finite temperature 
D^^^ ijth clement of D^^^ 

retarded propagator for spin 1 particles 
T forward scattering amplitude 

F^j, field tensor for electromagnetic field 
Fi^^ non-abehan field tensor for left (right) handed field Vi^r 
/be Bose distribution 
fpD Fermi distribution 
Qi^i, the metric tensor (1,-1,-1,-1) 

G° time-ordered free thermal propagator for nucleons 
Gp Hartree nucleon propagator; vacuum part 
Hartree nucleon propagator; medium part 

H^i, photon tensor 

J^(J^) hadronic (leptonic) electromagnetic current 

L^i, leptonic tensor 

Al invariant amplitude 

Mb Borel mass 

Mjv nucleon mass 

my mass of vector meson V 

M the invariant mass 

retarded improper self energy for spin 1 particles 

V vector mesons 

V three volume 

retarded electromagnetic current correlation function 
Z{13) partition function at temperature T — 1/(3 
a electromagnetic coupling constant, e = y/Ana 

as strong coupling constant, Qs = \/47ra7 

^aPui/ totally antisymmetric tensor, with e°^^^ = 1 
e thermodynamic energy density 

Fy width of the vector meson V 
ft four-volume 

retarded proper self energy for spin 1 particles 
Pixu spectral function for spin 1 particles 
Qfj^i, non-abelian field tensor for the p meson 

The subscripts L and T will be used to denote the longitudinal and transverse 
components respectively of a 2-ranked tensor in a medium. 
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1 Introduction 



The QCD renormalization group calculation predicts that strongly interacting sys- 
tems at very high temperature and/or density are composed of weakly interacting 
quarks and gluons [|I], 0, ^ due to asymptotic freedom and the Debye screening 
of colour charge. On the other hand at low temperature and density the quarks 
and gluons are confined within the hadrons. Therefore, a phase transition is ex- 
pected to take place at an intermediate value of temperature and/or density ||^, 
This transition is actually observed in lattice QCD numerical simulations |^ at 
high temperature. A system of thermalized strongly interacting matter where the 
properties of the system are governed by the quark and gluon degrees of freedom 
is called Quark Gluon Plasma (QGP). One expects that ultra-relativistic heavy 
ion collisions (URHIC) at CERN/SPS, BNL/RHIC and CERN/LHC might cre- 
ate conditions conducive for the formation and study of QGP |0, |^ (for recent 



development see Ref. UTOf). 

Various model calculations have been performed to look for observable signatures 
of this state of matter. However, among various signatures of QGP, photons and 
dileptons are known to be advantageous as these signals probe the entire volume 
of the plasma, with little interaction and thus, are better markers of the space- 
time history of the evolving fireball. This is primarily so because electromagnetic 
interaction is strong enough to lead to detectable signal, yet it is weak enough 
to allow the produced particles (real photons and dileptons) to escape the system 
without further interaction, carrying information of the fundamental constituents 
and their momentum distribution in the thermal bath. The real and virtual photon 
(dilepton) emission rate from QGP is determined by the fundamental theory of 
strong interaction, QCD. The dominant processes for the photon production from 
QGP are the annihilation {qq —* g'-f) and Compton processes {q{q) — > 9(^7)7)- The 
emission rate resulting from these reactions have been evaluated |jTl], [1^ in the 
framework of Hard Thermal Loops (HTL) resummation in QCD |T^. 



The disadvantage with photons is the substantial background from various pro- 
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cesses (thermal and non-thermal) [0, |T^, |T^, |T^. Among these, the contribution 
from hard QCD processes is well understood in the framework of perturbative QCD 
and the yield from hadronic decays e.g. 7r°/?7 ^77 can be estimated by invariant 
mass analysis. On the other hand, photons from the thermalized hadronic gas pose 
a more difficult task to disentangle. Therefore, to detect the electromagnetic signals 
from QGP it is very important to estimate photons (both real and virtual) from hot 
and dense hadronic gas along with the possible modification of the hadronic proper- 
ties below the critical temperature of the phase transition. However, the progress in 
our understanding of hot and dense hadronic matter has been retarded due to the 
nonperturbative QCD dynamics in the low energy regime. Because of this severe 
constraint considerable amount of work has been done on model building (see e.g. 



Refs. [^, ||2^, [^]) in order to study the low energy hadronic states. Never- 

theless, in URHIC, hadronic matter is expected to be formed after a phase transition 
from QGP. Even if such a phase transition does not occur, realization of hadronic 
matter at high temperature (~ 150 - 200 MeV) and/or baryon density (a few times 
normal nuclear matter density) is inevitable. As a result the study of hadronic in- 
teractions at high temperature and density assumes great significance. Also, there 
are several other aspects where medium effects may play an important role. For 
example spontaneously broken chiral symmetry of the normal hadronic world is ex- 
pected to be restored at high temperature and density and this will be reflected 
in the thermal modiflcation of the hadronic spectral function [^, |2^, ^ ^ 



These modiflcations can be studied by analyzing photon, dilepton as well as hadronic 
spectra. 

Various investigations have addressed the issue of temperature and density de- 
pendence of hadronic spectra within different models over the past several years. 
In particular, in-medium QCD sum rules are useful to make constraints on the 
hadronic spectral functions at flnite temperature and density p9|. Brown and 
Rho (BR) |]5D| argued that requiring chiral symmetry (in particular the QCD trace 
anomaly) yields an approximate scaling relation between various effective hadronic 
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masses, m*j^/mN ~ m*^/mfj ~ ''Tipi^/mp^^ ~ /t^/Ztt, which imphes, that all hadronic 
masses (except pseudoscalar mesons) decrease with temperature. In the Walecka 
model approach the vector meson mass gets shifted due to the decrease 

of the nucleon mass which appears through thermal loops in the vector meson self 



energy ||3^. The reduction in p meson mass has also been observed in the gauged 



linear sigma model @] at low temperature, however, near the chiral transition point 



it shows an upward trend. The nonlinear sigma model [^, |36[, claimed to be the 
closest low energy description of QCD shows the opposite trend, i.e. m* increasing 
with temperature. A similar qualitative behaviour of the p mass has been observed 



in the hidden local symmetry approach |]S^. The relation between the self energy 
and the forward scattering amplitude has also been utilized to study the change of 
hadronic properties in the medium ||38|, where the effects of non-zero temperature 



is rather small [3£]. Thus, there exists a lot of controversy in the literature about 
the finite temperature properties of hadrons. In view of this in the present article 
we shall consider various scenarios for the shift in the hadronic spectral function at 
finite temperature and evaluate its effects on the experimentally measurable quan- 
tities (as these issues should be settled by experiment), the photon and dilepton 
spectra originating from a thermalized system formed after URHIC. 

At this point we would like to clarify the scope of the present work. Our aim is 
to contrast the photon (both real and virtual) emission rate from the following two 
nuclear collision scenarios: 

(i) A +A ->QGP^Mixed Phase-^Hadronic Phase 

or 

(ii) A -|- A — ^Hadronic Phase, 

by taking into account the finite temperature effects on the hadronic masses and de- 
cay widths, calculated within the ambit of various models for hadronic interactions. 
The properties of matter formed after the nuclear collisions at SPS would most likely 
be determined by soft interactions, where the application of perturbative QCD is 
not comfortably justified and it is rather difficult to make any definite conclusion 
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about the formation of QGP here (for a recent review see [0). Therefore, at SPS 
we consider both the scenarios (i) and (ii). However, at RHIC the centre of mass 
energy of the colhsion increases by an order of magnitude where in addition to the 
soft processes the semi-hard processes contribute significantly to the initial energy 
density. We have considered scenario (i) exclusively for RHIC energies. This is be- 
cause, as shown later, the initial temperature estimated is too high for the hadronic 
degrees of freedom to survive. At even higher energies at LHC, the predictions of 
perturbative QCD become more reliable and the theoretical estimate for the ini- 



tial energy density is rather high |^ where the formation of QGP is beyond any 
reasonable doubt! 

The models considered to evaluate the medium effects are the Walecka model, 
the gauged linear and non-linear sigma models, the hidden local symmetry approach. 



BR scaling and Nambu scaling scenarios [42]. The QCD sum rules have been used 
to constrain the hadronic spectral function at non-zero temperature. Other models 



e.g, those proposed by Rapp et al and by Klingl et al |Q, where the effects of 
non-zero baryon density (baryonic chemical potential) is dominant over finite tem- 
perature effects are not discussed in the present article. This is because the aim 
of the present work is to study the hot, baryon free (central rapidity) region of the 
ultra-relativistic heavy ion collisions. As will be shown in the next section, both 
the photon and dilepton emission rates are proportional to the retarded electromag- 
netic current correlation function. These correlators or the spectral functions can 
be constructed in vacuum from the experimental data obtained in e^e~ hadrons 
(or from hadronic decays of the r lepton) for various isovector and isoscalar chan- 
nels 0. The in-medium spectral function of the vector meson is obtained, 
within the ambit of various models, by modifying the pole and the continuum struc- 
ture, resulting from its interaction with the constituents of the thermal bath. 

We organize the paper as follows. In section 2 we review the formalism of 
the emission of real and virtual photons from a thermalized system of strongly 
interacting particles. In section 3 we introduce the HTL resummation technique and 
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discuss the specific reactions which are used to calculate the photon and dilepton 
spectra. Section 4 is devoted to the discussion of the properties of hadrons at finite 
temperature. The gauged linear sigma model, the gauged non-linear sigma model 
and the hidden local symmetry approach have been described in section 5. The QCD 
sum rule approach has been discussed in section 6. We discuss space time evolution 
dynamics in section 7. In section 8 we present the results of our calculations and 
finally in section 9, we give a summary and outlook. 

2 Electromagnetic Probes - Formulation 

The importance of the electromagnetic probes for the study of thermodynamic state 
of the evolving matter was first proposed by Feinberg in 1976 While for most 
purposes one can calculate the emission rates in a classical framework, Feinberg 
showed that the emission rates can be related to the electromagnetic current-current 
correlation function in a thermalized system in a quantum picture and, more im- 
portantly, in a nonperturbative manner. Generally the production of a particle 
which interacts weakly with the constituents of the thermal bath (the constituents 
may interact strongly among themselves) can always be expressed in terms of the 
discontinuities or imaginary parts of the self energies of that particle ^|. In 



this section, therefore, we look at the connection between the electromagnetic emis- 
sion rates (real and virtual photons) and the photon spectral function ( which is 
connected with the discontinuities in the interacting propagators) in a thermal sys- 



tem ||50|, which in turn is connected to the hadronic electromagnetic current-current 



correlation function |5T| through Maxwell equations. It will be shown that the pho- 
ton emission rate can be obtained from the dilepton emission rate by appropriate 
modifications. 

We begin our discussion with the dilepton production rate. Following Wel- 
don let us define as the exact Heisenberg photon field which is the source 
of the leptonic current J' . To lowest order in the electromagnetic coupling, the 
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scattering matrix element, Shi, for the transition | /) H; l^l ) is given by 

Sni^e J {H;l+l- \ | /)rf^a;e*«^ (2.1) 

where | /) is the initial state corresponding to the two incoming nuclei, | H; 
is the final state which corresponds to a lepton pair plus anything, the parameter e 
is the renormalized charge and q = {qo,Q) is the four momentum of the lepton pair. 
Since we assume that the lepton pair does not interact with the emitting system, 
the matrix element can be factorized in the following way 

{H;l+l- I Jl{x)A^{x) I /) = {H I A'^ix) \ \ J^x) \ 0). (2.2) 

where | 0) is the vacuum state. Putting the explicit form of the current in terms 
of the Dirac spinors ('u(pi) and ^(^2)) we obtain 

Shi = e^^^^2^ [ d'xe--{H | A^{x) | /). (2.3) 



where 7^ denote the Dirac matrices, Ei — \Jpi +w?, with i — lh2 are the energy 
of the leptons and V is the volume of the system. 

The dilepton multiplicity N from a thermal system is obtained by summing 
over the final states and averaging over the initial states with a weight factor 

where Ei is the total energy in the initial state, Z[I3) is the partition function and 
(5 — 1/T is the inverse temperature. After some algebra N can be written in a 
compact form as follows: 

JS -e L ^M.^27r)3^i (27r)3£;2' ^^"^^ 
where L^^, is the leptonic tensor defined by 

spins 

= KP2+P^Pr-y5"^ (2.6) 



and H^y is the photon tensor 

H,, ^ -i- e-'^^« / d'xdSe^'-'^^-y^ E \ A,{x) AM I ^) e"^^-- (2.7) 

To obtain the above equation we have used the resolution of identity 1 = Yl,i I -^)(-^ I 
and the energy conservation equation Ei = + go, where go is the energy of the 
lepton pair and E^ is the energy of the rest of the system. Using translational 
invariance of the matrix element we can write 

H,, = ne-^^^^D>M, (2.8) 

where Q is the four volume of the system. Dj^^^ is the component iD"^]^ of the exact 
photon propagator D^,^ in the real time formalism of thermal field theory which has 
a (2 X 2) matrix structure (see appendix). The time ordered propagator is the (1, 1) 
component of D^u. In coordinate space it is defined as 

^Dllix) ^ ^^j:{H\T{A,{x)MO)}\H)e-^^" 

^ e{xo)D>,{x) + 9{-x,)D<M- (2.9) 

where D"^^ is the component iDj^l (see appendix), xo is the time component of 
the four vector, x^(= xq, x) and 9{xo) is the step function. Using the integral 
representation of the ^-functions in the above expression and taking the Fourier 



transform we get ||52| 



~ uj + ie (f — uj — ie 
Using the Kubo Martin Schwinger (KMS) relation in momentum space. 



(2.10) 



D>(a;,g-) = e^-D<,(^,g-), (2.11) 

we have 

^IM^. q) = -Y^^ImD;;,(go, q). (2.12) 
The rate of dilepton production per unit volume (N/Q) is then obtained as 

dN 2e^ d'pi d^P2 

d% = -^^^^ '"^^-^^"'^^M^ (2^- ^'-^'^ 
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Now, the spectral function of the (virtual) photon in the thermal bath is con- 
ventionally defined as 

P,u{q,,q) = / d'xe'^-Y.{H I [A,{x),AM] I H)e-f'^", (2.14) 



H 



SO that, we have ^ 



DfMo,Q)=l P^'^^^^^ -2i7TfBEiqo)pi,u{qo,q), (2.15) 
J—oo qo — u! -\- le 



where fBE^qo) = [^^'"^ — 1] ^- This leads to 

ImDlliqo, q) = -7r[l + 2fBE{qo)]pA%^ ^- (2-16) 
In terms of the photon spectral function the dilepton emission rate is obtained as 

O 2rf,u I ^ d?Pl d^P2 f I . 

— = 2vre p,.(go, q)j^^ J^J^^^^^^' (2.17) 

This relation which expresses the dilepton emission rate in terms of the spectral 
function of the photon in the medium is an important result. 

As is well known, it is not the time-ordered propagator that has the required 
analytic properties in a heat bath, but rather the retarded one. We thus introduce 
the retarded propagator which will enable us to express the dilepton rate in terms 
of the retarded photon self energy. The retarded photon propagator is defined as 

^DU1o^q) ^zW)l ^'^"^'■'^^(^o)?^^ I [M^)^MO)] I H)e-^^", (2.18) 
which leads to the relation 



H 



1. 



p,, = --ImD;,. (2.19) 



IT 

The above equation implies that in order to evaluate the spectral function at T 7^ 
we need to know the imaginary part of the retarded propagator. It is interesting to 
note that the above expression for spectral function reduces to its vacuum value as 
f3 ^ 00 since the only state which enters in the spectral function is the vacuum [p^]. 
Consequently we have 

1^ - -2e^L-I.«« (,..,-,^ (2.20) 
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This result can also be derived directly from Eq. ( p.l3| ) using the relation ^ 



lmDll = {l + 2fBE)lmD;^^. (2.21) 

Now the exact retarded photon propagator can be expressed in terms of the proper 
self energy through the Dyson-Schwinger equation: 

7-)R _ ^A"^ ^^f^^ _L _ 1 ^Zi^ ('9 99'! 

where, — «n^^ is the sum of all IPI (one particle irreducible) retarded photon self 
energy insertions which can be decomposed as 

nj, = vnf + 5^,nf . (2.23) 

Here A^^" and B^^'^ are transverse and longitudinal projection tensors respectively 
(see appendix) and 11^ and are the transverse and longitudinal components of 
the retarded photon self energy. The presence of the parameter ( indicates the gauge 
dependence of the propagator. Although the gauge dependence cancels out in the 
calculation of physical quantities, one should, however, be careful when extracting 
physical quantities from the propagator directly, especially in the non-abelian gauge 
theory. 

Inserting the imaginary part of the retarded photon propagator from Eq. ( p.22| ) 
in Eq. {^2^ we get 



|^ = 2.e=V(^-p. + B..,,)^^ (2.24) 

with 

1 imn^^i 

^^'"^ " ~^ (g2 + Ren^i)2 + (Imn^i)2- 



Comparing with Eq. ( |2.17| ) we have 

^ ^ (2.26) 

It has been argued by Weldon [55| that the electromagnetic plasma resonance oc- 



curring through the spectral function derived above could be a signal of the de- 
confinement phase transition provided the plasma life time is long enough for the 
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establishment of the resonance. Using the relation 



and 



2m^, / 4m2 
+ (2.27) 



q'^A^, = 

q'^B^, = 

g^'A,, = 2 

(^^'^i?^, = 1, (2.28) 



the dilepton rate is finally obtained as 



dR a n , 2rn? , / Am? , ^ „ , ^ , ^ 

A ^ "I^^^^ + " ^ PL)fBE{qo), (2.29) 

where m is the lepton mass and a denotes the fine structure constant. This is the 
exact expression for the dilepton emission rate from a thermal medium of interacting 
particles. In most of the cases the dilepton production rate from a thermal system 
is calculated with the approximation n^p = Ilf = 11^. Since ^lt both 

proportional to a, they are small for all practical purposes (this corresponds to the 
free propagation of the virtual photon in the thermal bath). Therefore pl^t{= p) is 
given by 

llmn^ llmnf'^ 
p = — = ^. 2.30 

^ IT IT 3q^ ^ ^ 



Using Eqs. ( ^1291) and (|O0|) we get 



dR a 2m^, L Am? ^ ^ , , 

= 12^7? ^ -rr - 

This is the familiar result most widely used for the dilepton emission rate It 
must be emphasized that this relation is valid only to O(e^) since it does not account 
for the possible reinteractions of the virtual photon on its way out of the bath and 
multiple emission of photon is ignored here. The possibility of emission of more 
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than one photon has also been neglected here. However, the expression is true to 
all orders in strong interaction. 

The emission rate of dilepton can also be obtained in terms of the electromag- 
netic current-current correlation function ||51|. Denoting the hadronic part of the 
electromagnetic current operator by J^, the leptonic part by and the free photon 
propagator by Dq'^ , we have the matrix element for this transition : 

SHi = e{H;l+r I j d'xdSJl{x)D^,\x-y)Jl:{y)\I). (2.32) 



This obviously follows from Eq. (|2.1|) by realizing that the solution of the equation 
of motion of the interacting photon field is 

A^^{x) = J dSD^'^ix - y)J%). (2.33) 

As in the earlier case the leptonic part of the current can be easily factored out. 
Writing the Fourier transform of photon propagator and squaring the matrix ele- 
ments, one obtains the rate of dilepton production 

^^-'^ ^>^^^{2^fE, (27r)3E,' (2.34) 
where VF^(g) is the Fourier transform of the thermal expectation value of the elec- 
tromagnetic current-current correlation function: 

W>M = j d'xe^^-j:iH I Jiix)j!:iO) I H) ——. (2.35) 



H 



The subtleties of the thermal averaging have been elucidated earlier. It is thus 
readily seen (Eq. ( p.34| )) that the dilepton (and photon) data yield considerable in- 
formation about the thermal state of the hadronic system; at least the full tensor 
structure of W^'^ can in principle be determined. Now, the current-current correla- 
tion function is related to the retarded correlator by 

W>, = 2e^^^fBEiqo)lmW;^u (2.36) 

where 

/p—PEjj 

dW-^(xo) E I [J»^ I ^) ^7^' (2-37) 



H 
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We then get 



We now define the (improper) photon self energy through the relation 



(2.38) 



(2.39) 



where —iP^y is the sum of all self energy diagrams. The advantage is that P/^ can 



R 



flU 



be defined in coordinate space as 



(2.40) 



Taking the Fourier transform and comparing with Eq. (p.37|) we see that P^iy{q) 
-iy^(g). Using the relations |]23|, and |]2| we end up with ^| 



dR 



a 



:i + ^)^^l^WlmP^' fBEiqo), 



q^ 



where 



q^ 



pRf, ^ gi^^pR = 2P^ + Pi 



(2.41) 



(2.42) 



It is important to realize that the analysis is essentially nonperturbative up to this 
point. To O(e^) we note that P reduces to the proper self energy 11 (= PDqD^^) 
and consequently Eq. (|2.41|) reduces to Eq. (|2.31|) . This approximation is the same 
as implied in Eq. (p.30| ). 

The connection between the electromagnetic current correlation function and 
the spectral function can be expressed in a straight forward way by substituting 
and using Maxwell equation {dad'' - ("^C - 1)^^* {da ^") = J^) in Eq. ( ^ ) 
to obtain 



2 C-1 2 C-1 

-q^qa hD> \q gpu - 



q 9m 



C 



C 

C-1 



q^qv j 

qfsq.] (l + fBE)- (2.43) 



Substituting the above equation in Eq. (|2.34|) we can recover Eq. ( p.lTj) . This 
establishes the connection between the two approaches of Refs.[|y, [5T[] . 
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The electromagnetic decay of unstable particles {e.g. p, uo and 0) within a ther- 
mal system could provide valuable information about the nature of the system. In a 
thermal medium the production of an off-shell vector meson iV) of four momentum 
q (where = M^) and its subsequent decay into a lepton pair leads to the dilepton 



emission rate as |^ 



dR = ^,pl.P'''fBE{qo)Tv^m- d\ (2.44) 

where p^^ is the spectral function of the off-shell vector meson given by Eq. ( p.l4| ) 
with the photon field replaced by the interpolating fields for vector mesons, the exact 
form of which is not important in the present case, P^u{= I] ^n^u) = —gfiu + Qfj.lu/l'^ 
is the projection operator for the vector meson V and Ty^i+i- is the partial decay 
width for the process V in vacuum. The spectral function is expressed in 

terms of the retarded vector meson propagator (as has been done before in Eq. ( p.l9| ) 
in case of photon). In the limit 11^ = 11^ = 11^, the spectral function is given by 

1 TmTT^ 



TT 



(g2 -ml + Relief + (Imn^)^ 



Using the relation P^'^ P^j^u = (2 J + 1), we get the dilepton emission rate due to the 
decay of an unstable vector meson of spin J as 

1 Imn^ 



(2J+1) 



(2.46) 



TT (g2 - + Ren^)2 + (ImU^y _ 

where Im 11^ is the imaginary part of the self energy of the particle V which should 
be calculated within the framework of thermal field theory ||59|. For a particle which 
does not decay in the collision volume (the total width Ftot = Imll^/ M is small) the 
spectral function in the above equation (term within the square bracket) becomes 
5(g^ — my), as it should be for a stable particle. In a medium the width Ftot for V 
should be calculated with all the processes involving the creation and annihilation 
of V, i.e. Ftot = Fy^aii - Faii^y [g. 

To obtain the real photon emission rate per unit volume [dR) from a system 
in thermal equilibrium we note that the dilepton emission rate differs from the 
photon emission rate in the following way. The factor -^^i// q'^ which is the product 
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of the electromagnetic vertex 7* — /+/^, the leptonic current involving Dirac 
spinors and the square of the photon propagator should be replaced by the factor 
J2 ^^i^t{= ~9fMu) for the real (on-shell) photon. Finally the phase space factor 
d^Pi/[{2nf El] d^p2/[{2TrfE2] should be replaced by d^q/[{2nfqo] to obtain 

dR = g^"" W> (2.47) 

2(27r)3^ ^"^ go 

As in the case of dileptons this expression can be reduced to 

This result can also be obtained directly from Eq. ( |2.31| ). The emission rate given 
above is correct up to order in electromagnetic interaction but exact, in principle, 
to all order in strong interaction. However, for all practical purposes one is able to 
evaluate up to a finite order of loop expansion. Now it is clear from the above results 
that to evaluate photon and dilepton emission rate from a thermal system we need to 
evaluate the imaginary part of the photon self energy. The Cutkosky rules at finite 
temperature or the thermal cutting rules 0, ^ give a systematic procedure to 
calculate the imaginary part of a Feynman diagram. The Cutkosky rule expresses 
the imaginary part of the n-loop amplitude in terms of physical amplitude of lower 
order (n — 1 loop or lower). This is shown schematically in Fig. (|1]). When the 
imaginary part of the self energy is calculated up to and including L order loops 
where L satisfies x + y < L + 1, then one obtains the photon emission rate for the 
reaction x particles y particles + 7 and the above formalism becomes equivalent 
to the relativistic kinetic theory formalism |^ . 



For a reaction 1 + 2^3+7 the photon (of energy E) emission rate is given 

by H 

= / ds / dt\M\^ / dEi 

d-^P 16(27r)'_E J(mi+m2)2 Jt^i„ J 

X / ,EJ-^^Mm±Im, (2.49) 

laE"^ + 26^2 + c 



where 



a = —{s + t~ml — mlY 
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Figure 1: Optical Theorem in Quantum Field Theory 
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M is the overall degeneracy of the particles 1 and 2, is the invariant amplitude of 
the reaction (summed over final states and averaged over initial states), / denotes 
the thermal distribution functions and s, u are the usual Mandelstam variables. 

In a similar way the dilepton emission rate for a reaction a a — > l~ can be 
obtained as 



dR 
'dFq 



d^Pa 



2Ea{27^y 



-J{Pa) 



d^Pa 



2Ea{2n) 



:f{Pa) 



d^pi 



d^P2 



2Ei{2tt)^ J 2E2(27r) 



M {2n)'S^'\pa+p-a-pi-p2)S^'\q-Pa-p-a). (2.50) 



where f{pa) is the appropriate occupation probability for bosons or fermions. The 
Pauli blocking of the lepton pair in the final state has been neglected in the above 
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equation. 



3 Photon and Dilepton Emission Rate 

In this section we briefly discuss the HTL resummation technique and the specific 
reactions considered in the present work to evaluate the electromagnetic probes from 
QGP as well as hadronic matter. 

3.1 Photon emission from QGP in the HTL approximation 

In the system formed after URHIC, numerous interactions take place among the 
charged constituents, either quarks or hadrons, which lead to the production of 
photons and dileptons. In this section we will discuss the emission rate of hard 
{E > T) photons from QGP. 

Naively, one expects that the properties of QGP at high temperature (T >> Tc) 
can be studied by applying perturbation theory due to the small value of the strong 
coupling constant, as{T). However, QCD perturbation theory at high temperature 
is plagued by infra-red problems and gauge dependence of the physical quantities, 
e.g. the gluon damping rate (see Refs. and ||66|). The gauge dependence 



of the gluon damping rate was cured by Braaten and Pisarski by an effective 
expansion in terms of hard thermal loops - i.e. including all the relevant loop 
effects in a given order of the coupling constant in a systematic way (for a beautiful 
review on gluon damping see Ref. [^). But the problem of infra-red divergences 
in QCD is not solved completely by the HTL framework. The quantities which are 
quadratically divergent in naive perturbation theory such as the damping rate of fast 
moving fermions in QGP becomes logarithmically divergent in effective perturbation 
theory. On the other hand quantities which are logarithmically divergent in the naive 
perturbation theory turns out to be finite if one applies HTL resummation method. 
The hard photon {E > T) emission rate which falls in the second category, is the 
relevant quantity for the present discussions. 

The thermal photon emission rate from QGP is governed by the following La- 
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grangian density: 



where 



^QGP — ^QCD + -^75, (3.1) 



^ a=l /=1 ^ 

1 

^ /=i 

In the above, is the non-abehan field tensor for the gluon field of color 
a, ipf is the Dirac field for the quark flavor /, Qg is the color charge, Cf is the 
(fractional) electric charge of quark flavor /, A"'s are the Gell-Mann matrices, F^y 
is the electromagnetic fleld tensor and A^^ is the photon fleld. As mentioned in 
the introduction the dominant processes for photon production from QGP are the 
annihilation [qq — » (77) and the Compton processes {q{q) However, 
the production rate from these processes diverges due to the exchange of massless 
particles. As mentioned earlier this is a well-known problem in thermal perturbative 
expansion of non-abelian gauge theory which suffers from infra-red divergences. 
One type of the divergences could be cured by taking into account the 'electric 



type' screening through the HTL approximation The non-abelian gauge theory 



also contains 'magnetic type' divergences, which can be eliminated if there is a 



screening of the magnetic fleld |6^, This is in sharp contrast to Quantum 
Electrodynamics, which is free from screening of static magnetic fleld. However, 
the study of magnetic screening is beyond the scope of HTL approximation as the 
transverse component of the gluon self energy vanishes in the static limit in this 
framework. Magnetic screening is relevant if any physical quantity is sensitive to 



the scale g^T, at which all the loop contributions are of the same order |70[ and 



hence the perturbation theory breaks down [71|. The production of soft photons 



{E < QsT) from QGP is non-perturbative because it is sensitive to the magnetic 



screening mass of the gluons []72| and consequently the soft photon emission rate is 
poorly known. Therefore, in the present work we consider only the production of 
hard photons {E > T) which is insensitive to the scale g^T. For such cases (hard 
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photon emission) the infra-red divergences could be ehminated within the framework 
of HTL as discussed below. 

The theory of HTL begins with the observation that at non-zero temperature 
there are two energy scales - one associated with the temperature T, referred to as 
the hard scale and the other connected with the fermionic mass ~ QgT [g^ « 1), 
induced by the temperature, known as the soft scale. A momentum appearing in 
the self energy diagram of photon would be called soft (hard) if both the temporal 
and the spatial components are ~ QsT (any component is ~ T). If any physical 
quantity is sensitive to the soft scale then HTL resummation becomes essential, i. e. 
in such cases the correlation function has to be expanded in terms of the effective 
vertices and propagators, where the effective quantities are the corresponding bare 
quantities plus the high temperature limit of one loop corrections. 

The notion of HTL can be clearly demonstrated in massless 0^ theory in the 
following way. Consider the Lagrangian density 

C = \{dct>f-g'ct>\ (3.3) 

The thermal mass (self energy) resulting from the one loop tadpole diagram in this 
model is m^j^ ~ Q^T"^- At soft momentum scale {p^ ~ gT) the inverse of the bare 
propagator goes as ~ g^T"^-- Thus, the one loop (tadpole) correction is as large as 
the tree amplitude. Therefore, this tadpole is a HTL by definition. Braaten and 
Pisarski |]14[ have argued that these HTL contributions should be taken into account 
consistently by re-ordering the perturbation series in terms of effective vertices and 
propagators. Therefore, according to their prescription we have 

C = ^{d<pr - 9^^' - \<<^^ + \<<^^ = ^eff + C^u (3.4) 

where Cct = '"^th*/'^/^ is the counter term which should be treated in the same footing 
as the 0*^ term. Cct has been introduced in order to avoid thermal corrections at 
higher order which has already been included in the tree level. With the counter term 
the Lagrangian remains unchanged, so the effective theory is a mere re-ordering of 
the perturbative expansion. A similar exercise has to be carried out in gauge theory 

21 




Figure 2: Two loop contribution to the photon self energy. A diagram interchanging 
the blob in the internal line of the third diagram should also be considered. 

keeping in mind that an addition and subtraction of local mass terms will violate 
gauge invariance. The effective action for hot gauge theories have been derived in 
Refs. 1?^, |73|, |7^, [76], |T^, whereas the authors of Refs. |T9[ follow the classical 



kinetic theory approach for the derivation of the HTL contributions. It has been 
shown in Ref. ||8^ that the contribution of HTL to the energy of the QGP is positive. 
The counter term required to avoid double counting in evaluating the virtual photon 
production from QGP in the two-loop approximation has been derived in Ref. [pi] 
recently. 

The photon emission from Compton and annihilation processes can be calculated 
from the imaginary parts of the first two diagrams in Fig. (H). Since these processes 
involve exchange of massless quarks in the t/u channels the rate becomes infrared 
divergent. One then obtains the hard contribution by introducing a lower cut-off 
to render the integrals finite. In doing so, some part of the phase space is left 
out and the rate becomes cut-off dependent. The photon rate from this (soft) 
part of the phase space is then handled using HTL resummation technique. The 
application of HTL to hard photon emission rate was first performed in Refs. [[ill , 
For hard photon emission, one of the (soft) quark propagators in the photon self 
energy diagram should be replaced by effective quark propagators (third diagram in 
Fig. d^)), which consists of the bare propagator and the high temperature limit of 
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Figure 3: Two loop photon diagram relevant for bremsstrahlung processes. The 
blob on the gluon (spiral line) indicates effective gluon propagator. The circle on 
the vertices represent those required to evaluate the imaginary part of the photon 
self energy in the framework of thermal cutting rules (see Refs. (|T^) and also ( [pT|] ) ) . 



one loop corrections When the hard and the soft contributions are added, 

the emission rate becomes finite because of the Landau damping of the exchanged 
quark in the thermal bath and the cut-off scale is canceled. The rate of hard photon 



emission is then obtained as 11 



= 5|^T^e--/-ln(2.912£Afr). (3.5) 

where is the strong coupling constant. Recently, the bremsstrahlung contribution 
to photon emission rate has been computed [|l3l by evaluating the photon self energy 



in two loop HTL approximation. The physical processes arising from two loop 
contribution (Fig. (^) are the bremsstrahlung of quarks, antiquarks and quark 
anti-quark annihilation with scattering in the thermal bath. The rate of photon 
production due to bremsstrahlung process for a two flavour thermal system with 
E > T is given by [|1^ 

= aasT' e"^/^ {Jt - Jl) In 2, (3.6) 

and the rate due to q — q annihilation with scattering in the thermal bath is given 
by, 

= ^ aa^ET e"^/^ ( - J,) , (3.7) 

where Jt ~ 4.45 and Jl ~ —4.26. The most important implication of this work is 
that the two loop contribution is of the same order of magnitude as those evaluated 
at one loop |n, |T^ due to the larger size of the available phase space. In case of soft 
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thermal photon {E ~ QsT) emission rate, all the vertices and the propagators has to 



be replaced by the corresponding effective quantities. It has been shown [0, ^ that 
the result is divergent due to the exchange of massless quarks introduced through 
the HTL effective vertices itself. However, such coUinear singularities for light-like 
external momentum could be removed with an improved action It is also shown 
that such infrared singularities could be removed (because of KLN (Kinoshita - Lee 



Nauenberg) theorem p6| , |8^ ) by including appropriate diagrams and summing 



over all degenerate initial and final states p8|, but the rate is non-perturbative 
because it is sensitive to the scale g^T ||72| . 

In this work, emission rate of hard photon is considered, which is well under 
control within HTL resummation. However, there are important issues in hot gauge 
theories which cannot be addressed within the HTL resummation method. For 
example, (i) HTL resummation is based on the weak coupling limit {qs << 1) to 
distinguish between hard (T) and soft momentum scale (gsT) but such a limit may 
not be realized in URHIC even for the highest energy to be available at the CERN 
LHC in the near future. Extrapolation of results obtained in HTL approximation 
to higher values of coupling constant will be demonstrated in section 8.2 through 
photon spectra, (ii) it cannot cure the infra-red divergence problem that arises in the 
damping rate of fast fermions, (iii)it cannot remove the mass shell singularities in 
the soft photon (real) emission rate, (iv) the next to leading order correction to the 
Debye mass diverges unless one includes magnetic screening, which is beyond the 
scope of HTL approximation and finally (v) HTL works for a system in equilibrium 
; extension of the formalism to non-equilibrium processes is still in the early stages 
of development. Results from other methods such as ladder approximation [pO| 



renormalization group equation etc. will be very important in these cases. For 
further discussions on the successes and limitations of HTL resummation technique 
and other methods we refer to Refs. |^ and 
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3.2 Photon emission from hot hadronic gas 

To evaluate the photon emission rate from a hadronic gas we model the system as 
consisting of vr, p, uj and rj. The relevant vertices for the reactions tttt — p7 and 
7T p TT'-f and the decay p n ii'-f are obtained from the following Lagrangian: 

C = -gpn.iT ■ (vf X d^n) - eJ^A^ + ^F^" {p^ x p,)^, (3.8) 

where Ffj^^ = d^Ay — dyA^, is the field tensor for electromagnetic field and is the 
hadronic part of the electromagnetic current given by 

= {Pu X ^^)3 + (7? X (9^7f + g,^^7: X /^))3, (3.9) 

with = d^pu — duPpL — gpmrip^i X Pu). vr, p and represent the tt, p and photon 
fields respectively and the arrows represent vectors in isospin space, gp^^^^ denotes 
the coupling strength of the p — n — tt vertex, fixed from the observed decay width 
p ^ 717T. The invariant amplitudes for all these reactions have been listed in the 
appendix of Ref. |Q . 

For the sake of completeness we have also considered the photon production due 
to the reactions nr] 7r7, tttt rj'-f and the decay u tt'-/ using the following 



interaction |92] 



£ = ^ e^ua^d^p^d-p^V + ^ e,.^f,d^uj''dy^n + ^A,p>^ (3.10) 
rriri gp 

where e^yap is the totally antisymmetric Levi-Civita tensor. The last term in 
the above Lagrangian is written down on the basis of Vector Meson Dominance 



(VMD) The invariant amplitudes for the reactions are given in Ref. jO^. The 
values of gppr, and g^p^, are fixed from the observed decays, p rj'y and ^ tt 7 
respectively ||9^. The constant gp is determined from the decay, p° e^e~ . 



The importance of the role of ai as an intermediary meson in the process vr p — >■ 
TT'-f was first emphasized in Refs. p3| , Recently it has been shown that the 



role of intermediary ai in this process is less important than thought earlier |9^, p^j . 



The photon production rate obtained in Ref. [p8| is similar to that in Ref. [p^ . 
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In this article we use the following interaction Lagrangian for the vrpai and nai'j 
vertices |8|, |9l |10( 



c 



z. 



m 



(2c + Z^)7r -p^x + ■ {d^p, - d,p^) x {d^a" - d'a^) 

■d^Ti ■ {d^pu - dyp^) X a" 



+ ^^^F^%d,a, - d^a, X 7f)3 (3.11) 

where corresponds to the ai field, is the renormalization constant for pion 
fields and fni= 93 MeV) is the pion decay constant. The interaction terms with 
coefficients c and are introduced to improve the phenomenology of the model. 
Following values of the parameters, = 1260 MeV, gp = 5.04, c = —0.12, Z^, = 
0.17 and = 1.25 |^9|, reproduce the width of the ai meson in vacuum. 



Most of the photon producing reactions under consideration involves unstable 
particles (e.g. p and uj) in the external lines (initial or final channels). According 
to Eq. (p.49|) we need to know the occupation probability of unstable particles in 
the thermal bath. So what is the appropriate phase space density of an unstable 
particle in this case? 

The density of a stable hadron of mass m in a thermal bath is completely de- 
termined by the temperature, chemical potential and the statistics obeyed by the 
species through the following equation: 

dN M 1 



d^xd^kds (2 vr)'^ exp(A;o — /i)/T ± 1 



6{s-m') (3.12) 



where J\f is the statistical degeneracy, kQ = \J + s is the energy of the particle in 
the rest frame of the thermal bath and p is the chemical potential. The question we 
would like to ask now - How the equation ( p.l2 ) will be modified if the particle decays 



within the thermal bath? This problem has been addressed by Weldon |^8|] through 
the generalization of Breit-Wigner formula at finite temperature and density. The 
distribution of an unstable particle in a thermal bath is given by 

dN M 1 



d^x d^k ds (2 exp(A;o - p)/T± 1 
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P{s) (3.13) 



where P{s) is the spectral function p.01|| and can be calculated from the effective 
thermal propagator. It is given by 

imn 



Pis) 



(3.14) 



n {s -ml + ReUy + (Imliy _ 
Equation ( p.l3| ) indicates that to obtain realistic results for the photon production 
through a reaction involving unstable particle in the external line the finite width of 
the particle should be taken into account by introducing the spectral representation 
of the corresponding particle and integrating over s. This is done in our calcula- 
tion for the unstable vector mesons appearing in the external line in reactions for 
photon production. In the case of an unstable particle appearing in the internal 
line the finite width of the particle is taken into account through effective propaga- 
tors. However, the effects of the finite width of p on the photon spectra is negligible 
but it affects the dilepton spectra substantially as will be shown in section 8. It 
is interesting to note that the spectral function reduces to a Dirac delta function 
6{s — rriy + Rell) in the limit Imll 0, i.e. when the particle is stable. In the 
calculation of the imaginary part of the self energy Imll of p say, one must in princi- 
ple, include all the processes which can create or annihilate a p in the thermal bath. 
However, within the ambit of the model adopted in the present work we have seen 
that the most dominant contribution to ImH comes from the p — 7i — n interaction 
in the temperature range of our interest. 

3.3 Dilepton emission from hot hadronic gas and QGP 

In order to express the dilepton emission rate from hadronic matter in terms of the 
retarded current correlation function we use Eqs. ( |2.27| ) and ( p.28| ) in Eq. ( |2.38| ) to 
obtain, 

where is the retarded current correlation function. The parametrized form 
of the electromagnetic current correlation function in the p and uj channels will 
be discussed in detail in section 5.2. Now instead of using the current correlation 
function directly in the above equation one can use vector meson dominance (VMD) 
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to obtain the dilepton yield from (tt+tt^ e+e^) which is known to be the most 
dominant source of dilepton production. In the low mass region one should also add 
the contributions from the decay of vector mesons such as p and u. This is usually 
done in the literature. In order to make a comparative study we state briefly how 
the emission rate from pion annihilation can be derived from Eq. (|3l5|) . VMD 
relates the hadronic electromagnetic current to the vector meson field through field 
current identity as 

= E (3-16) 

V 9v 

where, V stands for the vector fields p, u, 0. We shall keep only p meson in the 
following. The electromagnetic current correlator can then be expressed in terms of 
the propagator of the vector particle in the following way: 



1 4 

^ImD''^ 



(3.17) 



where 



111/ 



A 



flU 



(g2 - m2 + Ren^^)2 + [Imn^^]2 
ImUf 



(g2 -ml + Reli^l 



[Imnf]^ 



(3.18) 



Recall that A^^ and B^^ are the transverse and longitudinal projection operators 
whose explicit expressions are given in the appendix. In the approximation lif^ = 
we can define the p spectral function in VMD as 



ImW 



R 



2 4 

e ni p 
9l(f 



(g2 - m2 + Ren^-f^)2 + [\mI{pR\^ 



(3.19) 



Using Eq. (|3.15|) we obtain the dilepton rate from pion annihilation as 

2m2\ / 



dR 
dM 



M2 



Aw? 



imn^^ 



(M2 - m2 + ReUP^y + [Imn^«]2 
in the Boltzmann approximation. This on simplification gives 



dM (27r) 



n 



(3.20) 



(3.21) 
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where Ki is the modified Bessel function, and M is the invariant mass of the lepton 
pair and a{M) is the cross-section for the pion annihilation given by 

4 77 a' 



a{M) = — — ^1 - Aml/M^ ^1 - Am^/M^ (1 + 2myM^) \F^{M)\\ (3.22) 
3 



where 



4 

2 in 
In the same way, the invariant mass distribution of lepton pairs from the vector 



meson decays is obtained using Eq. ( |2.46|) , as 



^ = 2J + 1 K,{nM/T) 
dM 7r2 n n 

MFtot/vr 



-Mr™=_,,+ ,_, (3.24) 



tot 



where Ftot is the width of the vector meson in the medium and Fy^^+g- is the 
partial width for the leptonic decay mode for the off-shell vector mesons in vacuum 
given by 

TT-^e+e- = -^^Vl-4mVM2 (1 + 2mVM2) (3.25) 

where m is the mass of the electron. 

We have considered quark anti-quark annihilation for the evaluation of dilepton 
emission rate from QGP The dilepton rate is given by 

dM (27r)4 ^ n ^ ' 

with the cross section 



SOtt 



4 Hadronic Properties at Finite Temperature 

As emphasized earlier, the photon and dilepton emission rates are related to the 
imaginary part of the photon self energy in the medium. In this section we will study 
the in-medium modifications of the particles appearing in the internal thermal loop 
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of the photon self energy diagram. Here the hadronic medium consists of mesons 
and baryons at a finite temperature. Due to the interactions with real and virtual 
excitations, the properties of these hadrons are expected to get modified. As a result 
the propagators appearing in the photon self energy undergo modifications. The 
subject of the present section is to discuss how one incorporates these changes in the 
framework of Thermal Field Theory. A brief discussion of the thermal propagators 
is given in the appendix. 

4.1 The Walecka model - nucleon mass 

Before discussing the vector meson masses in the medium let us see how the nucleon 
properties are modified in matter at finite temperature. Nuclear matter is studied 



using the Quantum Hadrodynamics (QHD) model |^ in which the nucleons interact 
through the exchange of scalar a and the vector uj mesons. The interaction in this 
model is described by the Lagrangian 

Ci = -g^NNN-^^Nuj^" + g^NNNaN, (4.1) 

where N{x), a{x), and a;(x) are the nucleon, a, and u meson fields respectively. The 
cr{uj) field couples to the scalar (vector) current of the nucleon with the coupling 
constant QaNNiguiNN) which will be specified later. 

The free nucleon propagator at finite temperature and density in general has four 
components (see appendix). The time-ordered i.e. the (ll)-component is physically 
relevant for our purposes and we will denote this as G^{p). So we have 



N) 



+ 2m5{p^ - Ml)ri{p.u) 



p2 - Mjj + ie 

^ G'M + G'M, (4.2) 

where the first term {G%) describes the free propagation of nucleon-antinucleon pairs 
and the second term {G^^) allows for the on-shell propagation of particle-hole pairs. 
Mm in the above equation is the free nucleon mass. 
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Figure 4: Diagrammatic representation of Dyson- Schwinger equation for nucleons 
in RHA 



The effective mass of the nucleon in matter at finite temperature in the presence 
of interaction described by Eq. ([4.1| ) will appear as a pole of the effective nucleon 
propagator. In the Relativistic Hartree Approximation (RHA) |]3ll , 32] one obtains 
the effective propagator by summing up scalar and vector tadpole diagrams self- 
consistently i.e. by using the interacting propagators to determine the self energy. 
The effective propagator also called Hartree propagator is given by (see Fig. (^) 



G^(p)=G^(p) + G"(p)S^(p)G^(p) 



Hi 



(4.3) 



where S^(p) is the nucleon self energy which contains contributions from both scalar 



[T^s) and vector (S(^) tadpole diagrams [0, |3^ and is given by 



where 



and, 



yH _ ■ 



■ 'daNN 

1 

mi 



iJuiNN 



m 



Tr[G^(p)] 



2 / (2.)4Tr[7.G^(l>)] 



(4.4) 



(4.5) 



(4.6) 



Here, rrifj (m^) is the mass of the neutral scalar (vector) meson. The solution of 
Eq. ( ^I3| ) now reads, 

1 



G^(p) = (^+M 



N) 



p2 - ^ 



gUp) + g%{p) 



(4.7) 
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One observes that the pole structure of the effective nucleon propagator in RHA 
resembles that of the free propagator with shifted mass and four-momentum i.e. 
p = p + and = + T,^ , is the effective mass. Using in place of the full 
Hartree propagator in Eqs. ( |4.5| ) and ( [4.6| ) defines the Mean Field Theory (MFT) 
values of the self energies. This is equivalent to solving the meson field equations 
with the replacement of the meson field operators by their expectation values which 
become classical fields i.e. cr (cr) and u {u). This yields {a) = g^NN Ps/fn^ 
and {u'^) = Qujnn Pb I'mlj which indicate that the nuclear ground state contains 
scalar and vector meson condensates generated by baryon sources. The spatial part 
of the uo condensate vanishes due to rotational symmetry in infinite nuclear medium. 
These condensates are related to the scalar and vector self energies generated by 
summing tadpole diagrams in QHD as, = —gaNN{<7) and S° = —gu;NN{'^^)- 
The mean field approximation is thus to neglect the fluctuations in the meson fields 
which themselves are generated by the nucleons. 

RHA is obtained when one includes the vacuum fluctuation corrections to the 
MFT results. This amounts to the inclusion of the Dirac part of the propagator 
in the calculation of the self energies. Summing over the vacuum tadpoles 
results in a sum over all occupied states in the negative energy sea of nucleons. 
Vacuum or quantum fluctuations, as these are called, form an essential ingredient 
in a relativistic theory of many particle systems. Since there are infinite number of 
negative energy states in the vacuum one expects that the vacuum contribution to 
the self energy is infinite. 

Let us now find the Hartree self energy of the nucleon with the full nucleon 
propagator which consists a medium and a vacuum part. The vector part of the self 
energy is obtained from Eq. ( [4.6| ) as 



The first term of this equation appears to be divergent. The usual procedure is to 
regularize the integral in n dimensions by dimensional regularization to render the 
integral finite. One can then shift the integration variable from p to p. The resulting 




(4.8) 



32 



integral vanishes by symmetric integration. The vector self energy then reduces to 
= —gliNN^^^ PB/f^t thus gives rise to a an effective chemical potential, 
ji* = jj, — gtjNN Ps/f^lj- The scalar part of the self energy follows from Eq. 



where 



■9aNN 



mi 



(27r)4 ^2 _ + ie 



NN 



mi 



X 



fFD{p*,T) + fFD{p*,T) 



N 



(27r)3 E* 



(4.9) 



fFD{p\T) 

fFD{P*,T) 

E* 



exp[(E* + 1 

1 

exp[(E* +/i*)/T] + 1 



(p2 + m;2 



(4.10) 



Here Pb is the baryon density (in the present work we will take ps = 0) of the 
medium and is given by 

4 



PB 



(2vr) = 



d'p[fMP*,T)-fFDip*,T)]. 



(4.11) 



The first term in Eq. (^4.9| ), to be denoted by represents the contribution to 
the scalar self energy from the filled Dirac sea and is ultraviolet divergent. We will 
now proceed to renormalize this divergent contribution. The first step is to isolate 
the divergences through dimensional regularization. This gives 



^7^^^r(2-n/2)^^,3 



mi 



27r2 



-M 



N 



^r(2_n/2) + 3M^Ef + 3M^sf + Ef ) 

rZ 



(4.12) 



m^ 2tt^ 

since = + Sf^. The divergence in E^^-* now appears in the pole of the 
F-f unction for physical dimension n = 4. The counter terms needed to remove the 
divergent contributions from the loop corrections to the measurable amplitudes are 

4 



/:cT = E«nCr7n!. 

n=l 



(4.13) 



33 



Including the contributions from the counter terms the renormahzed self energy 
becomes 



where 



^CTC 



n=0 



.gaNN, 



(4.15) 



The coefficients (««) are fixed by defining a set of renormalization conditions. Since 
the scalar density (ipip) is not a conserved quantity the tadpole diagrams appear 
in the self energy. The tadpole contribution must vanish in normal vacuum (free 
space) i.e. {a)o = 0. This is ensured by the term aia in Cct- «2C"^ is the meson 
mass counter term which ensures that is the physical (measured) mass. Since 
the original Lagrangian of QHD |31| does not contain cr^ and cr^ terms, at the tree 
level, three and four point meson amplitudes must vanish. The last two counter 
terms in Eq. ( |4.13| ) are chosen to maintain this condition at zero external momenta 
for the a meson when nucleon loop corrections are included. We thus have 



(27r)' 



ttn = -i{-gaNNTin - 1)! 

Consequently the effective nucleon mass reads 



(4.16) 
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' (27r)3 
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Mn? 



- MUM*n - Mn) 
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(4.17) 



The solution of this equation gives the effective nucleon mass M^ as a function of 
temperature and baryon density. At zero baryon density it can be parametrized as 

8.941 



Mn 



M, 



N 



1 - 0.0264 



/ r(GeV) ^ 
V 0.16 , 



(4.18) 
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4.2 The vector meson mass 

In a medium meson properties get modified due to its coupling to nuclear excitations. 
This modification is contained in the meson self energy which appears in the Dyson- 
Schwinger equation for the effective propagator in the medium. The interaction 
vertices are provided by the Lagrangian 

CvNN = 9VNN (^Ni.T'^NV,^ - ^Na^^T-NdW^^ , (4.19) 

where = {cu'^,/?^}, is the nucleon field and Tq = {l,r}. r are the Pauli 
matrices, Ky is the nucleon-vector meson tensor coupling constant will be specified 
later. 

The lowest order contribution to the vector meson self energy is expressed in 
terms of the self-consistent nucleon propagator described in Eq. ( [4.7] ). This is given 
by 

U^'^ik) = -2zg'y^j, J ^Tr [V^ik) G^(p)r^(-A;) G^{p + k)] , (4.20) 

where represents the meson- nucleon vertex function obtained from Eq.( |4.19| ) 
and is given by 

r^(fc) = for uj 

T'^(k) = ^^ + i^a^-k^- for p (4.21) 

where cr'^" = |[7^, 7"]- The vector meson self energy can be written as a sum of two 
parts 

W{k) = Ii^;{k) + n^"(fc), (4.22) 



where 



Ii^;{k) = -2zgl.^^ I ^jT[T^ik)G^ip)T^i-k)G^ip + k) 

U^^^ik) = -2zgl.^^ J ^jT[T^{k)G^{p)T^{-k)G^{p + k) 
+T^ik) Gg(p)r(-A;) Gf{p + k) 

+ r'^(A;) G^{p)T''{-k) GUv + k)] . (4.23) 
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Up" is the vacuum polarization. This is a bihnear function of Gp and hence describes 
the correction to the meson propagators due to couphng to NN excitations. The 
NN pairs can be excited only if the four-momentum carried by the mesons is in the 
time-like region (/c^ > 0). Hence the shift in the mass of the vector mesons due to 
vacuum polarization is caused by processes like V — * NN V where N represents 
nucleons in the modified Dirac sea having an effective mass M^, smaller than what 
it would be in free space. We have seen that Up" causes a substantial negative shift 
in the masses of vector mesons. /^From Eq. (^4.23|) we have 

U^'^ik) - -2^0^ [ Tr[r^(;^ + Mj^)r-(;^ + ^ + Mj^)] 

IFtotci naive power counting it can be seen that this part of the self energy is ultravi- 
olet divergent and has to be renormalized. A few comments about renormalizability 
of the interaction given by Eq. (|4.19| ) is in order here. At very large momenta the 
propagator for massless boson ~ 0(A;~^), whereas for massive vector bosons it goes 
as ~ 0(1). This poses severe problems to the renormalizability of the theory with 
massive vector bosons. However, in a gauge theory with spontaneous symmetry 
breaking the vector gauge bosons acquire mass in such a way that the renormaliz- 
ability of the theory is always preserved. The theory which involves neutral massive 
vector bosons coupled to a conserved current is also renormalizable. This is because 
in a physical process the propagator Dq'^ = {—g^" + p'^p'^ /m^)/ {p^ — ra^ + ie) ap- 
pears between two conserved currents and J^, and the offending term p^p^ jrn? 
does not contribute because of current conservation {p^J^ = 0), making the theory 
renormalizable. This is the case for the u meson which we shall consider first (see 
Refs. [rU^, |103|| ). The counter term required in this case is 



^y^iv = -^C^"^V;.- (4.25) 

We use dimensional regularization to separate the divergent and the finite parts. The 
divergences now appear as a pole in the gamma function at the physical dimension 
n = 4. The renormalized vacuum polarization tensor for the u is then given by 

U'P'ik) = [g^"" - PA;7fc2)H5,^"(fc2), (4.26) 
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where 



UujNN 

7r2 



r(2-n/2) / dzz{l-z) 



- dzz{l - z) ln[M;;,2 _ ^2^^i _ ^)]| _ ^ (4 27) 



in which the counter term contribution 



n 



liuCTC 



(4.28) 



has been included. C, is now determined by the renormahzation condition 



nr(fc')l 



0. 



(4.29) 



Finally, we arrive at 



Re(nr)i: 



^'^^-A;2 f' dzz(l-z)\n 







- k^zil - z) 



Ml - k^z{l - z) 



(4.30) 



Renormahzation of the vacuum self energy for the p meson presents additional 
problems because of the tensor interaction. A phenomenological subtraction proce- 



dure, as described in Refs. ||104| , |105|| is used to arrive at the following expressions: 



where 



_9pNN_p 



2M, 



Jo 



N 



dz z(l — z) In 



2 '2M, 



N 



M*^^ - k^ z(l - z) 



Ml - k^ z(l - z) 



dz In 



k'^z(l-z) 



Ml 



(4.31) 



(4.32) 



(4.33) 



k'^z{l-z) 

The medium dependent part of the polarization, II^J^, describes the coupling 
of the vector mesons to particle-hole excitations. It contains at least one on-shell 
nucleon propagator which provides a natural ultraviolet cutoff in the loop momenta. 
This part of the self energy leads to an increased effective mass of the vector mesons 
in the medium. 
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The change in the hadronic mass in the medium can be understood from the 



following phenomenological arguments ||106|| . Consider the propagation of a vector 
meson in a nuclear medium. The attenuation of the amplitude at a distance z, in 
a Fermi gas approximation, is given by e""'^^, where n is the density of nucleons 
and a is the meson-nucleon interaction cross section. The optical theorem relates 
(J to the imaginary part of the forward scattering amplitude; a = 4TTlmJ-'{E)/k. 
It then follows that the meson wave function ip ~ exp[27TinzJ^{E)/k]. In terms 
of an effective mass {rries = m + Am), the propagation can also be described by 
ip ~ expli^jE"^ — mlffz]. Comparing the arguments of the exponential we get 



Am 



m 



-ReJ^(E). 



(4.34) 



This relation clearly shows that the enhancement or reduction of hadronic masses 
depends on the sign of ReJ-'{E). 

In a hot and dense medium because of Lorentz invariance and current conserva- 
tion the general structure of the polarization tensor takes the form 

W^-" = Uriko, IkDA'^" + IlLiko, 1^1)5^^ (4.35) 

where the two Lorentz invariant functions 11^ and 11^ are obtained by contraction: 



n. 



(4.36) 



is the four velocity if the thermal bath. 



In the case of the vector meson interacting with real particle-hole excitations in 
the nuclear medium these are given by 
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-2tg 



VNN 



Tr 



(2vr)4 



T''{k)GF{p)T''{-k)GD{p + k) + {F^D) 



(4.37) 



with 



(n 



9vNN 1 f Pdp 
27r2 i Up 

- 8\p\\k 



(r + 2M*^') In 



{k' + 2\p\\k\?-<u;'p 
(P-2|p1|fc|)2-4A;2a;2 



(4.38) 
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UJr, 



(A;2-2|p1|A;|)2-4A;2cu2 



X 



(4.39) 



9vNN (j^y ^ f 

47r2 V2Mjv/ J 



k f pdp 
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(k^ - 2\p\\k\)^ - Ak^o<^^ 
fpDil^*,T)+fMl^*,T) 



The longitudinal component of the polarization tensor is given by 



(4.40) 
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(4.42) 



(4.43) 



(4.44) 
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In the above the superscripts 'f', 'ft' and 't' represent the vector- vector, vector- 
tensor and tensor-tensor components respectively arising from the product of vector 



and tensor terms in Eq. ( [4.21|) . The dispersion relation for the longitudinal (trans- 
verse) mode now reads 

kl-\k\'^-ml + Renf(y)(fco, k) + ReU^{k'^) = (4.45) 

Usually the physical mass (my) is defined as the lowest zero of the above equation 
in the limit k 0. In this limit 11;^ = = U^, and we have, 



In^ = n = + 

3 ^ 



(4.46) 



where 

U^{ko,k-^0) -- 
with 

F{\P\,M*^) 
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VNN 



p'dpFm,M*^)[fpn{f^*,T) + fFD{fi\T)] (4.47) 
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lim' + 3M*^') + ki {2Mu^: 



N 



(4.48) 



where Up=p^ + M'ff. 



The effective mass of the vector meson is then obtained by solving the equation: 



kl-ml + Ren = 0. 



(4.49) 



The effective masses (denoted by asterix) take the following parametrized forms: 

5.24" 



rrin 



1-0.127 



1 - 0.0438 



T(GeV) 

Y(GeV) 
0.16 



7.09' 



(4.50) 



The effective mass of the ai meson, m*^ has been estimated from m* by using 
Weinberg's sum rule ||107| . 
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One finds reference to two otlier kinds of masses in the literature. Tlie invariant 
mass is defined as tlie lowest order zero of Eq. ( [4.45|) with 11^ neglected. Again, the 



screening mass of a vector meson is obtained from the pure imaginary zero of the 
quantity on the left hand side of the same equation with ko = 0. These two masses 
are different because of the non-analyticity of the polarization tensor at the origin 
i.e. at {ko,k) = (0,0). 

5 Models with Chiral Symmetry 

In this section we will discuss finite temperature effects on the vector meson prop- 
erties for those models which respect chiral symmetry. The effects of in-medium 
properties of vector mesons on the electromagnetic ejectiles will be presented in 
section 8. 

5.1 The gauged linear sigma model 

The linear sigma model (LSM) is a beautiful tool to describe the low energy dy- 
namics of pions, because it shows explicitly how the spontaneous symmetry breaking 
(SSB) of global chiral symmetry (SU(2){g)SU(2)) by the isosinglet a field generates 
pions, the Nambu-Goldstone (NG) bosons. However, there are reservations about 
the description of the cr meson as a well-defined degree of freedom because of its 
large decay width which is comparable to its mass. But, it has been argued by 
Hatsuda and Kunihiro (see also Ref. ||108| |) that in the limit of chiral symmetry 



restoration, the decay of a to two pion state should be disallowed as cr and tt become 
degenerate in mass in this limit. They have explicitly shown that the width of a due 
to (T ^ 27r decay vanishes as T ^ T^, where is the critical temperature for chiral 
transition (it is still not known whether the critical temperature for deconfinement 
and chiral transition are the same or not, in the present work no distinction is made 
between them). 

The properties of vector mesons at finite temperature within the ambit of gauged 
LSM are studied by Pisarski |25, 3i]. In the following we will discuss the main results 
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of his work. 

The simplest version of LSM contains isosinglet a field and isotriplet pion field 
and respects the charge conjugation, parity and time reversal symmetry (CPT) [p5|] . 
The Lagrangian obeying these constraints is 

CLSM = tr I 5^$ |2 +^2^r I $ |2 +]-\tr{\ $ \''f - htri^) (5.1) 



where $ is defined as 



$ = ^(a + ivf ■ f) (5.2) 



with r being the Pauli matrices. The non-zero value of h ensures that the pions 
are massive and consequently PCAC (partially conserved axial current) relation is 
satisfied. Note that | $ p= + tt^ is chirally invariant and elimination of the a 
field by imposing the condition cr^ + tt^ = results in the Non-Linear Sigma Model, 
which will be discussed in the next section. 

In the gauged LSM (see Refs. and ||109|| for detail discussions) one introduces 
the vectors and their chiral partners (axial vector) through left and right handed 
fields as follows, 

= ifT + an ■ t + {uj'^ + fn (5.3) 

\//^ = (p^ - ■ t + (cu^ - /f ) (5.4) 

where Oi and /i are the chiral partners of the p and u mesons respectively and 
t = t/2. 

The inclusion of vector (axial) mesons would increase the number of possible 
couplings and hence the number of arbitrary parameters become large. How to 
include vector meson in LSM with minimal coupling to the matter fields (vr and 
cr)? To perform this, following Kroll et al ||11(]|| , one assumes that the Lagrangian 
and its chiral transformation properties should be such that the current generated 
by the chiral transformation is proportional to the vector field itself, which leads 
to the relations known as field current identities. This requirement automatically 



leads to the idea of VMD due to Sakurai |93]. The field current identity is achieved 



by promoting the SU(2) global chiral symmetry of vector fields to a local gauge 
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symmetry as was done for the first time by Yang and Mills ||1 1 1|| for the isospin 
symmetry. The Lagrangian for the vector field reads, 

A. = \tr I Fr ? +\tr I Fr T +niltr[{Vrf + (K^)^] (5.5) 

where Fj^!^ = Vi^r~^^ ^i^r~^9 ^!r; • above Lagrangian the kinetic term 

for the gauge fields remains invariant under the transformation and the field-current 
identity is obtained through Gell-Mann Levy theorem from the mass term of the 
gauge fields as 

Jtr = "f^l'r (5.6) 

We note that chiral symmetry is a global one in QCD, therefore, the local symmetry 
has to be broken and this is achieved precisely by the mass term of the vector fields in 
the Lagrangian. Next one has to introduce the interaction of matter fields (vr and a) 
with the gauge fields preserving the field current identity. Noting that the ordinary 
derivatives occurring in Eq. ( p.l|) spoils the field-current identity, we introduce the 
required interactions consistent with the gauge principle i.e. by replacing the partial 
derivatives by covariant derivatives: 

Di'^ = d^'^-ig{Vl'^-^VI') (5.7) 

Finally the resulting Lagrangian density for the gauged LSM is obtained from 
Eqs. (13), (U) and as 

Cgism = tr I D^^ \^+iiHr I $ P+^Atr(| $ - htr{^) + Cir (5.8) 

Expanding the kinetic term for the matter field one finds 
1 



tr I 



|2_ _ 



{d,a + ga ■ + {d,ir + gp^ x tt - ga'^a)' + g\a' + Ti')U^f 

(5.9) 

The above equation indicates that (i) a shift in the a field (o" — ^ ctq + o") gives rise 
to mixing between vr and ai fields (a term ~ ga^d^Tt ■ arises from the second 
term of the r.h.s. of Eq. ( ^.9|) ), which has to be eliminated by an appropriate shift 
in the ai field, (ii)there is no interaction term involving which can come into the 
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picture through anomaly and (iii) the kinetic term for pion gets modified because of 
the shift in the ai field. Thus to get back the cannonical form of this term one has 
to renormalize the pion field vr n/^/Z^, where Z^, = rn?p/rn?^^] Z-J^ = 1/2 gives 
the Kawarabayashi - Suzuki - Riazuddin- Fayyazuddin (KSRF) relation ||112| , |113 



After some algebra one gets ml = h/{Z.„ao), ml = h/ao + 2AcrQ, = y/Z^ao. 
Taking m^ = 137 MeV, m„ = 600 MeV, = 770 MeV and = 1260 MeV, we 
get do = 152 MeV, h = (102 MeVf, /i = 412 MeV and A = 7.6. 

With these inputs Pisarski has evaluated the thermal masses of p and ai at 



lowest order in g at low temperature. We quote the results below p4 



Ami (S'^p + 4A;2) 



im 



m 



2\2 



Km 



2 

Tfl A- 

45m2 



P L 



Ami (3m^ + Ak"^ 

(3,1 V Cti ' 



2ml 



m 



ai 



(5.10) 



(5.11) 



im^ — m^l^ ml [ml —ml) ml 
In the chiral limit (ctq goes to zero and many of the couplings vanish) assuming 
the validity of VMD in the medium Pisarski has showed that ||114| p and ai become 



degenerate with a mass value ~ 962 MeV (p mass increases). On the other hand, 
if one adopts a scenario where vector meson dominance (VMD) is not valid in the 
medium then mp(T^) = ma^iT^) = 630 MeV (p mass decreases). However, it is 
important to mention at this point that the chiral symmetry can also be realized 



via the Georgi limit [115] where the p meson becomes massless. Pisarski |25] has 
argued that the results obtained by Georgi in the non-linear sigma model can be 
translated in terms of gauged linear sigma model without the validity of VMD, 
for which there is no unique prediction for the behaviour of p mass at non-zero 
temperature. Thus the behaviour of in-medium p depends on the validity of VMD 
in the medium. The thermal shift of the pole of the spectral function of p in the 
gauged LSM are shown in Fig. (0) of section 8. 

5.2 The gauged non- linear sigma model 

It is well-known that the global SU{2)i SU{2)r symmetry of 2 fiavours QCD is 
expected to be spontaneously broken to the subgroup SU{2)v and the pions appear 
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as the N-G boson. The non-hnear sigma model with SU{2)i^ SU{2)r / SU{2)v is 
an effective theory of QCD for the description of pion dynamics. The in-medium 
properties of vector mesons have been studied by Song |116|| in the framework 



of gauged non-hnear sigma model (NLSM) [^. This model will be discussed only 
briefly because it is very similar to the gauged LSM: the main difference is that the 
a degree of freedom is eliminated in NLSM by the non-linear realization of chiral 
symmetry as mentioned in the previous sections. We start with the observation that 
a perfectly valid parametrization of $ could be, 

u = ^M^i:%]^-M^<P] (5-12) 

where = 0"r"/A/2 is the pseudoscalar field and = ^J^f^^. The Lagrangian for 
the NLSM based on the manifold 5t/(2); SU{2\ISU{2)y is given by, 

J IT . 



Co = ^tr 



d^Ud^U^ (5.13) 



The vector and the axial vector fields can be introduced as the Yang-Mills gauge 
as before to minimize the number of arbitrary parameters in the model. The result- 
ing Lagrangian is given by. 



Cnlsm — 



-tr I Fr I' --tr I Fr \^ 



+mltr[{Vn' + {Vn'] (5.14) 



where V"^'^ = {v^ ± a^)/2, and denote vector and axial vector fields. To 
improve the phenomenology of the model, the following higher dimensional terms 



can be added to the Lagrangian [^, |117|| without spoiling the symmetry under 
consideration, 

jO-mm = -t^tr [d^UDM^ F''^'' + D^U^ F^'^"] (5.15) 

where ^ is a constant determined form the decay of vector mesons [Q. The thermal 
shift of the p-mass has been evaluated by Song with pion loop, pion tadpole and pion- 
ai loop resulting from the interaction terms of the Lagrangian given in Eqs. ( ^.14| ) 
and ( [5.15[ ). It shows negligible change in the p-mass from its vacuum value. 
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The effective masses of p, ai and u at non-zero temperature liave also been 
calculated by Song ||116|| with a SU (3)/ SU{3)r symmetric Lagrangian: 



NLSM 



4 



D^uDf'u^] - ^tr I rr \^ P 



+m^ir[(\/,'^)2 + (V;'^)^] + -/^r [M(f/ + f/t - 2) 



(5.16) 



where ?7 is defined as U = exp[— ^'^ i 



2i 



exp[ 



and Aa's are Gell-Mann 



matrices. The two higher dimensional terms with co-efficients ^ and k are added 
to improve the phenomenology. Please note that although these terms retain the 
gauge invariance of the model, they spoil the renormalizability of the model. 

The dynamics of uj is governed by the anomalous interaction, also known as 
Wess-Zumino interaction given by 



anomaly 



(5.17) 



This is very similar to the Gell-Mann Sharp Wagner interaction already encoun- 
tered in section 3.2. 



Song ||116|| has considered the following values of the parameters consistent 
with the vacuum properties of the vector and axial vector mesons: {g,n,C,) = 
(10.30,0.34,0.45) and (6.45, -0.29,0.06), referred to as set I and II respectively. 

The calculation of thermal mass shift of the vector and axial vector mesons with 
all these inputs reveal that: (i)with parameter set I p and u masses increase with 
different rate and ai mass decreases, (ii) with parameter set II the thermal mass shift 
of p and u is negligibly small but ai mass decreases slightly. The modification of 
hadronic masses due to thermal interactions within the ambit of the model discussed 
above are presented in Fig. (|I^) of section 8 . 

5.3 The hidden local symmetry approach 

In case of the two chiral models described above the vector mesons are introduced as 
Yang-Mills field and the mass term for the gauge boson are put in by hand which may 
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not be entirely satisfactory, whereas in the hidden local symmetry (HLS) approach 
the p meson is generated as a dynamical gauge boson of a hidden symmetry in the 
NLSM ||118| , |1 19|| . It has been explicitly shown by Bando et al that, in general. 



any NLSM corresponding to the manifold G/H is gauge equivalent to a "linear" 
model having Ggiobai Hiocai symmetry. Accordingly, the Lagrangian of Eq. ( ^.131 ) 
can be written in a form that exhibits, besides SU{2)i^ SU{2)r global, a local 
SU{2)v symmetry - the hidden symmetry and p meson appears as a gauge boson 
corresponding to this symmetry. (The axial vector ai is not included in the minimal 
version of HLS Lagrangian.) To make it more explicit, we introduce two SU{2)- 
matrix valued variables C,i{x) and C,r{x) with the transformation properties |118|| , 

^i,r{x) h{x)^i,r{x)gl (5.18) 

with 

U = ij^r (5.19) 
where h{x) E [SU{2)v]iocai and gi^r e [S U (2) i^r] global- 6,r is parametrized as 

6,r = exp[zS(x)//sT (5.20) 

where vr = vr'^t" and S = S^t*^. We note at this point that the unwanted degrees of 
freedom, S, which have entered in the system via Eqs.( |5.19| ) and ( p.2(]| ) are known 



as "compensator"- the would be N-G boson which has to be eaten up by the hidden 
gauge boson, p and these extra degrees of freedom reappear as the longitudinal 
polarization of the (massive) p. Now we define the covariant derivative as 

D,^i = d^^i-tgV,^i + t^il, (5.21) 

and 

D^^r = d^ir - igVf^^r + i^r T (5.22) 

where l^{r^) is the external field corresponding to the gauging of SU{2)i<^ SU{2)r 
and is the gauge field corresponding to the symmetry [SU{2)y]iocai- With 
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these fields two types of [SU{2)i^ SU{2)r]giobai®[SU{2)v]iocai invariants can be 
constructed |lT|, |12|, 

J n 



Cy = -ftr 



D,^i ■ ^1 + D^^r ■ 6] (5.23) 
D,^i-^}-D^^r-^lT (5.24) 



and a linear combination C = Ca + ctjCy is equivalent to the original Lagrangian 
given in Eq. ( |5.13| ). By fixing the gauge ^/ = = exp(i7r//7r) (and hence eliminating 
the unphysical degrees of freedom, S) one can show that Ca = Cq, while Cy vanishes 
when equation of motion for is used. So far has been treated as an auxiliary 
field. It is assumed that the kinetic term for this field is generated by quantum 
effects or by the QCD dynamics. The full Lagrangian with the kinetic term is 

jChls = Ca + aCv - ^g^,u0^'' (5.25) 

where Qf^i, is the non-abelian field tensor for the p meson. The Lagrangian of 
Eq. (|5.25|) can be written as ||118|| , 



Chls = lid.n)' + ^agfT ■ n x d,n + ^g'afy^ - ^g,,^'' + (5.26) 



The above equation implies that the mass of the p meson (m^ = ag'^f^) is generated 
due to SSB via Higgs mechanism and the unphysical N-G modes (S not vr) are 
"eaten-up" by the gauge boson i. e. the three extra degrees of freedom get converted 
to the three longitudinal polarization of the massive gauge boson. For a = 2 one 
recovers the KSRF relation. This value of a also results in universal coupling of p. 



Harada et al [37] have evaluated the finite temperature effects on the p-mass 
upto one loop order in the HLS approach due to the thermal pion and p meson 
interactions. Their results reveal that at high temperature the reduction in p mass 
due to pion loop is overwhelmed by the increase due to thermal p loop contribution, 
although the net shift is rather small. The contribution of thermal pions to the p self 
energy in this model is different from other calculations because in HLS approach 
there is no pion tadpole contribution. We will study the effects of the mass shift of 
p in this approach to the electromagnetic probes later in this article. 
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6 Spectral Constraints at Finite T 



In the previous section we have discussed the change in the hadronic properties 
using effective Lagrangian approach. However, medium modifications can also be 
studied by applying QCD sum rules (QSR) [ |121| , |122| , |123| |. As many good reviews 



are available on the QSR at zero temperature |46, |124| , 125 , 126 , 127 after the 
original work of Shifman et al ||128|| , we, therefore, very briefly introduce the basic 
principles of QSR in vacuum and then discuss the QSR at non-zero temperature. 

6.1 QCD sum rules at zero temperature 

The basic aim of the QCD sum rule approach is to evaluate the resonance parameter 
(mass, coupling constant etc) in low energy hadronic physics in terms of the vacuum 
expectation values of quantities such as quark and gluon condensates, (0|gg|0) (in- 
dicates chiral symmetry breaking), (0|G'^^G'^'^"|0) (signals the break down of scale 
invariance), etc, where g(x) is the quark field and G^^j^i^x) is the gluon field tensor. 
These condensates are non-perturbative in nature and appear as a power correc- 
tions to the leading logarithmic (perturbative) behaviour. These power corrections 
are far more important than higher order corrections [p,24||. In the following we 



will discuss how the QCD sum rule approach connects the perturbative and non- 
perturbative domain and leads to the determination of resonance parameter {i.e. 
mass and the resonance strength) of the p meson. 

The QCD sum rule approach starts with the Wilson operator product expan- 
sion for the time ordered product of two (or more) currents. The gluon and quark 
condensates appear as higher dimensional operators in the expansion. The coeffi- 
cients of this expansion contain the short distance part and the long distance part is 
contained in the vacuum expectation values. The coefficient can be evaluated per- 
turbatively in terms of the parameters (a^ and the quark masses) of the Lagrangian 
used. We consider the time-ordered or causal current correlator 

n^.(g) = I J rfVe*^-^(0|T(^J^(x)J,(0)) |0) 

= (g^^. - gV)n(g') (6.1) 
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In this article we will confine our discussions for currents with quantum numbers 

J', = \ («7m« - dl,d) (6.2) 

where u and d are the quark fields for up and down quark respectively. The analytic 
structure of the correlator (11), for spacelike = — g^, can be expressed through a 
dispersion relation: 

^/^9n 1 f Imn(s) (is , , , , , 

n g2 = _ / + subtraction 6.3 

TT J S + Q'' 

Imaginary part of 11 is proportional to the spectral density which can be modeled 
as consisting of a conspicuous resonance and a continuum with a sharp threshold 

Imn(s) =nY.gn ml 5{s - m|) + ^ f 1 + ^) 0(s - ^o) (6.4) 

Res 8^ ^ ^ ^ 

with a resonance strength Qji and a pole position at m^. 

The theoretical side of the sum rule is derived from an operator product ex- 
pansion for large = —q^ (deep Euclidean region where asymptotic freedom is 



realized) as has been suggested by Shifman et. al ||128|| . Thus we write 



I J d^xe"'-^T(^J,{x)Ux)) = C,(g) + E Cn(g)a (6.5) 

where / is the identity operator, C's are the Wilson coefficients, and 0„'s are the 
local gauge invariant operators constructed from the quark and gluon fields. The 
operators are ordered by their increasing dimensions and therefore, the coefficients 
fall off by the corresponding power of q^. On dimensional ground one sees that the 
operators of dimension d > Q leads to 1/q'^ power corrections. However, for large 
= —q^ a fewer power corrections {d = 6) is sufficient to converge the series. 
Taking the vacuum expectation value of Eq. (|6.5| ) we obtain [|124|1 



where 



n^,,. = (q^qu - q fi'/.i.)n (6.6) 



1 / a \ 1 
n = -—{l + — )\n^ + -—{0\m^uu + mddd\0) 



/i2 2Q4 
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1 rv TTfy — 



TTCy — 

^jO\(uj^X-u + dj,X^d) Yl ll.X^m (6.7) 

^ q=u,d. 



The left hand side (Lh.s.) is the well known two point Greens function which 
can be expressed in terms of the phenomenological parameters, characterizing the 
strong interaction processes, consistent with the current under consideration, via the 
dispersion relation. The right hand side (r.h.s.) has been evaluated by using OPE 
in the short distance (asymptotic freedom) region. The vacuum expectation value 
of the higher dimensional operator appears as a power correction to the asymptotic 
contribution (the first logarithmic term in the r.h.s. of the above equation). 

The sum rule therefore, becomes (modulo subtractions) 

1 I lmUis)ds ^ ^ 
n J s + 

In Eq. (|6.8| ) r.h.s. corresponds to large or small distance scale with fewer 
power corrections and l.h.s should be saturated by the lowest resonance, which is 
a long distance phenomenon. Therefore, in order to get a balance between the two 
sides we would like to have a weight function which enhances the low contribution 
relative to the high one. This can be done by taking additional derivative with 
respect to Q^, and then taking and the number of derivatives n to infinity, 
we obtain Borel transformed sum rule. Borel transformation is equivalent to the 
following mathematical operation: 



1 1 



where 



e-'/^'" (6.9) 



Lm = Urn Q'^ --— (6.10) 

Q2/^ — jV/^— const. 

and Mb is the Borel mass. Applying Eq. (|6.9| ) on the l.h.s. of Eq. (|0|) and Eq. ( |6.10| ) 
on the r.h.s. and expressing vacuum expectation value of four fermion operators in 
terms of two fermions, we obtain [ p. 24 

^2 



j e-''^B\raT[{s)ds = 



l + - + ^(OKgg|0) 
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+ 



448 TT^a, 2' 
^ mF^ '''' ^ 



(6.11) 



Substituting the various values of the matrix elements as given in Ref. ||124 
obtain 



we 



J e-'I^^B Imn(s) ds = ^M] 



a. 0-04 



0.03 



Differentiating Eq. (|6.12|) with respect to 1/M^ to obtain another sum rule: 



1 + — 

TT 



0.04 0.06 

+ 



M% 



(6.12) 



(6.13) 



In Eqs. |6.12| and |6.13| the terms and arise due to gluon and quark conden- 
sates respectively. Assuming that r.h.s. of Eq. ( |6.4D is saturated by the p resonance 
we get from Eqs. (|6.12|) and ( |6.13|) 



Ml 



(1 + a./vr) [1 - (1 + cJo/M|) e-^^l^^l \ - 0.04/M| + 0.06/M6 



(6 14) 

(1 + a,/7r) [1 - e-'^o/*'^s] + 0.04/M^ - 0.03/M| ^ ' ' 

The above expression still depends on the Borel mass Mb and the continuum thresh- 
old loq. The value of loq can be inferred from the data of e^e~ annihilation. The 
absolute value of p mass is then obtained by looking for the stability plateau i.e. 
choosing such that dmp{Ml) / dMi = 0. To determine the resonance strength 
for the p meson we keep only the p resonance in the sum of Eq. ( |6.4|) and substitute 
it in Eq. ( |6.12| ) to obtain, after an elementary integration, 

0.03 



27rm^ 



, as 0.04 

1 + — H r 

M% 



(1 + a,/7r)e-^«/^^- 



(6.15) 



Qp is related to Qp as Qp = 1/g't 



Eqs. ( |6.14| ) and ( |6.15|) indicate how the resonance parameters of vector mesons 
can be extracted by using QCD sum rules in vacuum. In the next section we will 
briefly discuss the QCD sum rules at non-zero temperature. 

6.2 QCD sum rule at non-zero temperature 

As mentioned earlier, the retarded correlator has the required analytic properties in 
a thermal system. QCD sum rules for vector mesons in the medium ||122| , |123|| start 
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with the retarded current correlation function (defined as W/^,, in section 2) 

Il^Mo,q)=tJ d'xe'^''e{xo){[Mx),JM]) , (6-16) 

where = {qo,Q) is the four momentum, with the source (electromagnetic) currents 
defined in terms of the quark fields as (in units of e) , 

2_ 1 - 1_ 

Jfi = -^u-f^u - -d-f^d - -s-f^s (6.17) 

Defining the current in the p, uj and channels as 

j; = (l/2)(«7^«-(i7^d), (6.18) 

r^^{l/2){u-f^u + d-f^d), (6.19) 

and 

= S7^s, (6.20) 
one can express the electromagnetic current in terms of p, u) and fields as, 

^. = ^.^ + ^'/;-^'/^ (6.21) 

As discussed earlier there are two independent invariants in the medium; the trans- 
verse (11;^) and the longitudinal (Hf ) components of the polarization tensor. In the 
limit g ^ 0, as there is no spatial direction, 11;^ and becomes equal (= n-^) and 
the trace of the retarded correlation function can be expressed in terms of as 
= nj[J^/(— S^o)- Both the transverse and the longitudinal components satisfy 
the fixed q dispersion relation. In particular, at g = 0, 

Ren^(go) = -P r du"^^^^-^ + (subtraction). (6.22) 

Rell'^ can be calculated using perturbation theory with power corrections (operator 
product expansion (OPE)) in the deep Euclidean region Qq — > — oo. For exam- 
ple, OPE for Ren-'^(go), which is the same as the OPE for the causal (Feynman) 
correlator 11^ (go), has a general form at = — — > — oo, 

Ren^(Q^ ^ -oo) ^ -Colng^ + g C^asin^Hf^VQ')) ^^^^^.^^^ ^^ ^S) 



n=l 
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where /i is the renormahzation point of the local composite operators, If \Q\ is 
much larger than other soft scales present in the system such as Aqcd and T, 
perturbative QCD can be used to calculate the relevant quantities. C„ are the c- 
number Wilson coefficients which are T independent. All the medium effects are 
contained in the thermal average of the local operators C„. Since {On)j. ~ T^'-Agji,^ 
with I + m = n due to dimensional reasons, (|6.23|) is a valid asymptotic expansion 
as long as ^ and Aq^j^. The local operators in the vector meson 

sum rule are essentially the same with those in the lepton-nucleon deep inelastic 
scattering (DIS) and can be characterized by their canonical dimension (d) and the 
twist (r=dimension-spin). They are given in Ref. | |123[ up to dimension 6 operators 



and we will not recapitulate them here. For q —>■ 0, Eq. (|6.23| ) is an asymptotic series 
in l/Q"^ or equivalently an expansion with respect to d. The medium condensates 
may be evaluated by low energy theorems, the parton distribution of 
hadrons and lattice QCD simulations. 

Matching the l.h.s. and the r.h.s. of Eq. ( |6.22|) in the asymptotic region 
gp — > — oo is the essential part of QSR. This procedure gives constraints on the 
spectral integral and hence the hadronic properties in the medium as well as in the 
vacuum. There are two major procedures for this matching, namely the Borel sum 
rules (BSR) [ |128|| (discussed in the previous section) and the finite energy sum rules 
(FESR) ||129|| , which can be summarized as 

I dql W{ql) [Imn«(go)-ImnJ,,(go)] = 0, (6.24) 

W(n^\ _ f Qo"" 0(^0 - qi) (FESR), 
^Wo) - j e-^oA^I (BSR). 

Here Imn^^^(go) is a hypothetical imaginary part of 11^ obtained from OPE and 

Mb is the Borel mass. 

We have seen earlier that in QSR in the vacuum, the spectral function (i.e. Imll^ 

in Eq. ( |6.22| )) is modeled with a resonance pole and the continuum to extract the 

mass and decay constant of hadrons. In the medium, such a simple parametrization 

is not always justified because of the thermal broadening of the spectrum and also 

because of the new spectral structure due to Landau damping and the thermal 
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mixing among mesons. Therefore, the model independent constraints obtained from 
QSR are only for the weighted spectral integral. 

For example, the first three finite energy sum rules at finite T read ||123 

i"00 

h = I [lmU\qo)-lmU^^^{qo)]dql = 0, (6.25) 
h = I [Imn^(go)-ImnJ^^(go)]go'^9o = -C^2(02),, (6.26) 



oo 



h = [\mn%q,)-\mnf{q,)]qtdq',=C,{0,),. (6.27) 







Similar sum rules hold for the axial vector channel (in the chiral limit) except that 
one has a different operator for O^^. One can also generalize the above sum rules to 
finite (fpOj [131|]. 



Explicit forms of Cn{On)j, have been calculated as [|123 



1 (y 

Co = --(1 + ^), Ci=0, (6.28) 

C^m, = ^{^G'), + Usq^loDoq),, (6.29) 

16 

^3(^3)7 = -(scalar 4 - quark) J, + —(5gz7o-Do-Do-Dog)r- (6.30) 

Here we have neglected the terms proportional to the light quark masses (chiral 
limit) and the quark-gluon mixed operators. The operator S is used to make the 
operators symmetric and traceless. At low T, one may use the soft pion theorems 
and the parton distribution of the pion to estimate the r.h.s. of the above equations. 
When T is close to T^, one has to look for a totally different way of estimation: the 
simplest approach is to assume the resonance gas to evaluate the r.h.s., while the 
direct lattice simulations will be the most reliable way in the future. An important 
feature of the OPE in the above is the appearance of local operators with Lorentz 
indices. This happens because we are taking the rest frame of the heat bath which 
breaks covariance. 

The sum rules can be used to check the validity of the calculations of the 
spectral functions using effective theories of QCD. This is in fact quite useful for the 
spectral function at finite baryon density. At finite T, especially near the critical 
point, the behaviour of the condensates with dimension d > 4 is not known precisely. 
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Therefore, it is rather difficult to make a strong argument on the spectral constraints 
near at present. The future lattice simulations of these condensates are highly 
called for. 



6.3 Parametrization of the spectral functions 

In this section we will introduce a parametrization of the correlator at finite T. 
The parametrization should be consistent with the experimental data from e+e" 
hadrons processes at zero T, and it should be also consistent with the high energy 
behaviour known from perturbative QCD at go ^ T. 

As the vector mesons appear as resonances in the electromagnetic correlator, 
using Eqs. (|6.16| ) and (|6.21|) we can write. 



ImnJ, = Imn^',« + ^Imn-;,^ + hmU%^. (6.31) 

The above equation shows that the contributions of u and mesons to the electro- 
magnetic probes are down by almost an order of magnitude compared to p meson. 

As mentioned before, at zero three momentum the imaginary part of the trace 
of the retarded correlator can be written in terms of its longitudinal component as, 

Im nJ,(go) = -3 Im nf (go) = -3go' Im flf (go). (6.32) 

Thus our next task is to parametrize Imllf (go). We take a Breit-Wigner form with 
an energy-dependent width for the resonance along with a continuum: 

At zero T, this reduces to a relativistic generalization of the parametrization used 
by Shuryak to ffi the experimental data of e^e^ hadrons. Here Dp is the 
imaginary part of the self-energy which should in principle contain all the channels 
which can destroy or create a p in the thermal bath. Hence Dp is given by the dif- 
ference of the decay-width and the formation width so that Dp = gor(go). However, 
we have seen that for a baryon free matter the most dominant contribution to Dp 
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comes from the pion-loop [ |132|| . For a p meson propagating with energy go and three 
momentum q the p width is given by 



rp-^7r7r(Q'0) Q) 



487r ^ ' go 



1 + 



2T 



X In 



l-exp[-f(go + Vr(s)Jg2-s 



l-eM-Uqo-W{s)Jql-s)] 



(6.34) 



where s = — (f and W{s) = yl — 4m^/s. In the hmit |g| — > , the above 
expression reduces to the in-medium decay width and is given by 



r,^_(go) = ^goiy^(go) 



1 + / 



BE\ 



1 + /. 



.go- 



BE\ 



fBE{^)fBE{'^) 



(6.35) 

with fBE^x) = [e^ — ^^0 is the continuum threshold above which the asymptotic 
freedom is restored and fp is the couphng between electromagnetic current and the 
p field defined as 

(0 I I p) = fprn^e, (6.36) 
Assuming vector dominance in the medium we obtain, 



9p = mp/ fp 



(6.37) 



In the vacuum, the standard parameters for the p spectral function are given by, 
rup = 0.77 GeV, = 0.14 GeV, fp = 0.141 GeV, gp = 5.46, coq = 1.3 GeV, S = 0.2 
GeV and as = 0.3. The resulting spectral function for the p-meson in the vacuum 
should be compared with Ref. P:6[ |. 

Let us now concentrate on the spectral function in the u channel. We again take 
a Breit-Wigner form along with a continuum: 

D,. 1 a.o, 1 



Im nf'"(go,g-'=0) = /^ 



where f^j is a coupling of the current with the cj-meson defined as 



(6.38) 



{0\ Jp \uj) = f^m^ep. 



(6.39) 
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Note that /^^ here is defined as factor 3 larger than Shuryak's definition . D, 
which is the imaginary part of the self-energy, is calculated using the Lagrangian 
density given in Eq. ( p.lOp . We have shown in earlier calculations ||132| , |133|| that 



a substantial contribution to the uj width comes from the process oon ^ tttt in a 
thermal bath. Consequently 

= 9o(r<^-*37r + r^^TT^TTTr) (6.40) 

where 

r.^3.(go) = C dw dx\ F \^ S (6.41) 

S is the phase space factor for thermal equilibrium, given by 

^ = [(1 + fBEm)ii + fBEm){i + fsEm) - fsEmfBEmfBEm] 

(6.42) 

and 

C = ^"^"^r"o^° (6.43) 
487r3 ml 

The limits of integration are 



^min 


= m^, 


^max 


= {q',-3ml)/2qo, 


•^max 


= -y/0.5u; (W - W7max 


a'lnin 


a^max; 


El 


= W, 


E2 


= x + {qo-w)/2, 


E3 


= -x + {qo-w)/2, 



\p^\ = \/Ef-ml, (6.44) 
and Pi is the pion 3-momentum. The amplitude for the process is 

|Fp=|p, iMl-^o)^ (6.45) 



where 



uj^ + ml- 2uj{Ei + E2) + 2E1E2 

Zo = — ^ T, ^ — 6.46 

2 I P1P2 I 
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and 



with 



H = j2K (6.47) 



=1 



h2 



h. 



3 



1 


?i2 + ml 




1 




Qi3 + ml 




1 




Qh + ml 





hi = 2(gi2gi3 + m^r^)/ii/i2 

^5 = 2(gi3523 + "2jrJ)/l2/l3 

K = 2(gfi2g23 + mlTl)hihz 

gi2 = (^i + ^2)'-p1-m2 

g23 = (E2 + Es)^-f,-ml (6.48) 

The width for a;7r ^ tttt is calculated analogously. 

In the vacuum the standard parameters for cu are as follows: m^^ = 0.782 GeV, 
= 0.14 GeV, = 0.138 GeV, coq = 1.1 GeV, S = 0.2 GeV and a, = 0.3. 
In a medium at finite T, we simply replace rUp, ujqJp and Qp by the corresponding 
effective quantities (denoted by asterix) m*, uJq, f* and g^* respectively. Since not 
much is known about the critical behaviour of the scalar and tensor condensates 
at finite T in QCD sum rules we take a simple ansatz for in-medium quantities for 
their T-dependence. A possible parametrization of *-quantities at finite T is 

= (6.49) 
mv fv i^o V ^cj 

where A is a sort of dynamical critical exponent and V stands for vector mesons 
{p and uj). (Note that there is no definite reason to believe that all the in- medium 
dynamical quantities are dictated by a single exponent A. This is a simplest possible 
ansatz.) Since the numerical value of A is not known, we take two typical cases: 
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A = 1/6 (BR scaling) and 1/2 (Nambu scaling) |^2|. The effective mass of ai is 
estimated by using Weinberg's sum rules [|107|| . 

Some remarks are in order here: 
(i) Eq. ( |6.49| ) for m* is not entirely consistent with the low temperature theo- 
134], which says there should be no O(T^) correction to the mass. There- 
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fore, one cannot take the ansatz too seriously at low T. In practical applications, 
however, T < 100 MeV is not relevant in any way since it is below the freeze-out 
temperature. 

(ii) Local duality constraint Ji in QCD sum rules implies that (/p)^ = 87r^(l + 
as/7f){ujQ)^ + (scattering term) [ p.23| ]. This condition is slightly violated for /* in 
Eq. (|6.49|) because of the existence of the scattering term (Landau damping). 

(iii) The assumption of vector dominance in the medium together with Eq. ( |6.49| ) 
simply leads to g*p = Qp- 

Under these reservations, we will use the parametrized spectral functions (BR 
scaling and Nambu scaling) in the calculation of the lepton and photon productions 
in later sections. Major qualitative difference between the spectral function in the 
effective Lagrangian approaches and that in this section is the existence of the 
continuum and its medium modification at finite T. 

As mentioned in the introduction the photon and dilepton emission is deter- 
mined by the retarded correlator of the electromagnetic current. Dilepton emission 
involving the p and uj mesons is thus obtained by inserting times Eqs. (|6.33|) and 
( |08|) in Eq. (g35|) using Eq. ([O^ ). 



7 Evolution Dynamics 

As mentioned earlier, in URHIC the produced matter will either be in the form of a 
hot hadronic gas or a quark gluon plasma. So far we have talked about the rate of 
photon and dilepton emission per unit time from unit volume of a thermal system 
made up of quark matter and hadronic matter at a fixed temperature T. Our next 
task is to consider its evolution in space and time. This is done using relativistic 
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hydrodynamics. A basic ingredient of the hydrodynamic description of the colhsion 
volume is the existence of a strong interaction time scale, 

Ti ~ -r-~ ~ Ifm/C ~ Tformation (7.1) 

In any hadronic collision the produced fragments can only interact after a proper 
time Ti has elapsed after their collisions. Thus, there is another time scale in the 
problem, the so called transit time , which is defined as 

Ttransit ~ (7.2) 

Ra is the nuclear radius, 7cm is the Lorentz factor. If the value of 7cm (which 
is a function of the collision energy) is such that Ttransit < ^formation then most of 
the secondaries are formed after the nuclei pass through each other. Consequently 
these secondaries will not contribute to the energy density of the fluid in the central 
region. Such a scenario may be realized at RHIC (Relativistic Heavy Ion Collider) 
and LHC (Large Hadron Collider) energies. This particular feature has been taken 
into account in Bjorken's hydrodynamic model ||135 . 



7.1 Bjorken's hydrodynamical model 

It has been observed experimentally that the particle spectra for the secondaries 
produced in — collisions exhibit a central plateau in the rapidity space. This 
kind of behaviour is due to the frame independence symmetry of the hydrodynamic 
expansion of the system ||136|| . Bjorken assumed that the same kind of plateau will 



also be observed in nucleus nucleus collisions ||135|| . In terms of the initial condition 
this means that the energy density, pressure etc (all the thermodynamic quantities) 
will be a function of the initial thermalization (proper) time only and will not 
depend on the space time rapidity t] (defined later). This initial symmetry of the 
thermodynamic quantities is preserved throughout the evolution scenario. If the 
particle rapidity density is flat or invariant under Lorentz boosts then the entropy 
density (s) will be independent of the rapidity. Since our discussion is limited to 
the baryon free region, there is only one independent thermodynamic variable T, 
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say. Once s is independent of Lorentz boost (rapidity) so are all the thermodynamic 
quantities. 

The evolution of the fluid is governed by the energy momentum conservation 
equation 

d^Tf"' = (7.3) 

where T^"^ = (e + P)u'^u'^ + g'^'^P is the energy momentum tensor for ideal fluid. 
For an isentropic flow the entropy conservation reads 

df,s^ = (7.4) 

where s'^ = su^ is the entropy current. Let us consider the frame independence 
symmetry in a two dimensional sub-space [t ~ z plane). Changing the independent 
variables from (t, z) to (r, r/) using 

1 , t + z 



= 7^23^. (7.5) 



the equation of motion Eqs. ([7.3|) and ([7.4|) become 



d d 

— {sT cosh(?/ — 7])) + ^ ("^ sinh(?/ — 77)) = (7.6) 

— {Tt smh{y -r])) + — {T cosh(y - r/)) = (7.7) 

The independent variable r, by definition is the proper time of the frame which is 
related to the cm. frame by a Lorentz transformation along the z-axis with velocity 
z/t. The variable t], known as the space time rapidity, becomes equal to the fluid 
rapidity y{= \ ln(l + v^) /(I — f z)). Putting y = r] in Eqs. ( |7.6| ) and ( [r.7| ) we get 

|- (st) = (7.8) 
dT 

These equations imply that T is independent of rj and so are all the thermodynamic 
quantities and st =const. This is the Bjorken's scaling solution. The resulting 
space-time picture of the collision is shown in Fig. (^. 
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Figure 5: Space-time diagram of the collision in Bjorken hydrodynamics 



It may be noted that the above results were obtained without any specific input 
from the equation of state (EOS), it is thus a general result that one dimensional 
similarity fiow is necessarily isentropic even if there is a phase transition. For a 
relativistic massless gas with statistical degeneracy g^, s and T are related through 
the equation of state: 

s = ^^9kT' (7.10) 

Putting this expression for entropy density in the Bjorken scaling solution we get 
T^T =const. This is the cooling law which is extensively used to evaluate the 
signals of QGP. The initial temperature of the system is determined by observing 
that the variation of temperature from its initial value Tj to final value Tf (freeze- 
out temperature) with proper time (r) is governed by the entropy conservation 
(Eq. m) 

s{T)t = s{T,)t, (7.11) 

The entropy density is then expressed in terms of the observed particle (pion) mul- 
tiplicity. Using Eqs. ( [7.10[ ) and ( |7.11| ) one gets the initial temperature as 

T.3 = ^ (7 12) 

where dNT^/dy is the total pion multiplicity, Ra is the radius of the system, Xj is the 
initial thermalization time, and C(3) is the Reimann zeta function, ak = gk 
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is the degeneracy of the produced system and hence k stands for either QGP or a 
hot hadronic gas. The rapidity density for the secondaries is obtained from ||137| , 



where a is known as the rescattering parameter. dN/dy |p_p can be parametrized 
to fit the experimental data in the central region as a function of the centre of mass 
energy, as 

(f) =°-«'"(^' P-"' 

\ ^ / p—p 

The assumption of a central plateau in the rapidity distribution is not experimentally 
observed in nucleus nucleus collisions at the presently available energies. Hence the 
boost invariant hydrodynamics may not be a valid concept at these energies. The 



concept of complete stopping in Landau model ||138|] is not valid either at these 
energies. The physical situation may be in between the boost invariant model of 
Bjorken and the Landau model of complete stopping, which means that there may 
be an overlap between the formation zone and the collision zone. The Bjorken 
model will be used in this work to describe the space time evolution of matter 
formed in URHIC. Appropriate generalization has been made to take into account 
the temperature dependent hadronic masses. 

7.2 Initial conditions and equation of state 

The set of hydrodynamic equations is not closed by itself; the number of unknown 
variables exceeds the number of equations by one. One thus needs to postulate a 
functional relation between any two variables so that the system becomes determin- 
istic. The most natural course is to look for such a relation between the pressure P 
and the energy density e, as is done in the case of thermal equilibrium. Under the 
assumption of local thermal equilibrium, this functional relation between P and e is 
the Equation Of States (EOS). Obviously, different EOS's will govern the hydrody- 
namic flow quite differently ||139|| and as far as the search for QGP is concerned, the 
goal is to look for distinctions in the observables due to the different EOS's (corre- 
sponding to the novel state of QGP vis-a-vis that for the usual hadronic matter). 
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It is thus imperative to understand in what respects the two EOS's differ and how 
they affect the evolution in space and time. Recently, the sensitivity of the photon 
emission rate on various evolution scenarios has been studied in Ref. ||140| . 



A physically intuitive way of understanding the role of the EOS in governing 
the hydrodynamic flow lies in the fact that the velocity of sound = {dP/de)^ 
sets an intrinsic scale in the hydrodynamic evolution. One can thus write a simple 
parametric form for the EOS: P = Cg(T)e. Inclusion of interactions, however, may 
drastically alter the value of cj. [|141|| . In our calculation we assume the MIT bag 
model equation of state for the QGP where the energy density and pressure are 
given by 

^Q=9Q^ + B, (7.15) 



and 



PQ = gQ^T'-B. (7.16) 



The effective degrees of freedom in QGP, qq = 37 for two flavours. The entropy 
density sq is given by sq = 2gQ{n^/4:5)T^. Putting ak = = {7i'^/90)gQ the initial 
temperature for a system produced as QGP can be determined from Eq. ( |7.12| ). 

In the hadronic phase we have to be more careful about the presence of heavier 
particles and the change in their masses due to finite temperature effects. The 
ideal limit of treating the hot hadronic matter as a gas of pions originated from 
the expectation that in the framework of local thermalization the system would 
be dominated by the lowest mass hadrons while the higher mass resonances would 
be Boltzmann suppressed. Indirect justification of this assumption comes from the 
experimental observation in high energy collisions that most of the secondaries are 
pions. Nevertheless, the temperature of the system is higher than during a major 
part of the evolution and at these temperatures the suppression of the higher mass 
resonances may not be complete. It may therefore be more realistic to include higher 
mass resonances in the hadronic sector, their relative abundances being governed 
by the condition of (assumed) thermodynamic equilibrium. We assume that the 
hadronic phase consists of vr, p, u, t], ai mesons and nucleons. The nucleons and 
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heavier mesons are expected to play an important role in the EOS in a scenario 
where mass of the hadrons decreases with temperature. 

The energy density and pressure for such a system of mesons and nucleons are 
given by 

^H= E / d'pEhfBE{Eh,T) + I d'pENfFD{EN,T) (7.17) 

h=mcsons V / V / 

and 

n=mesons ^ ' " ^ ' 

where the sum is over all the mesons under consideration and N stands for nucleons 



and Eh = y + ml ■ The entropy density is then 

SH = '-^^^ ^ 4aeff(T)T^ = A^g,sim*iT),T)T' (7.19) 

where g^s is the effective statistical degeneracy. Thus, we can visualize the finite 
mass of the hadrons having an effective degeneracy gcs{'m*{T),T). Because of the 
temperature dependence of the effective degeneracy Eq. ( |7.12| ) has to be solved 
self consistently in order to calculate the initial temperature of the system initially 
produced as a hot hadronic gas. We thus solve the equation 

dN^ 45C(3) 2 3 

= RA^^oE{Ti)TiTi (7.20) 

where aefr(Tj) = (tt^/QO) 5'efr('^*(^i)? • The change in the expansion dynamics 
as well as the value of the initial temperature due to medium effects enters the 
calculation of the photon emission rate through the effective statistical degeneracy. 

If the energy/entropy density in the fireball immediately after the so-called "for- 
mation time" Tj is sufficiently high, then the matter exists in the form of a QGP. 
As the hydrodynamic expansion starts, the system begins to cool until the critical 
temperature is reached at a time tq. At this instant, the phase transition to the 
hadronic matter starts. Assuming that the phase transition is a first order one, the 
released latent heat maintains the temperature of the system at the critical temper- 
ature Tc, even though the system continues to expand; the cooling due to expansion 
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is compensated by the latent heat hberated during the process. Together with the 
possible explosive events, we are neglecting the scenarios of supercooling or super- 
heating. This process continues until all the matter has converted to the hadronic 
phase at a time r//, the temperature remaining constant at T = T^, from then on, 
the system continues to expand, governed by the EOS of the hot hadronic matter 
till the freeze-out temperature Tf at the proper time tj. Thus the appearance of 
the so called mixed phase at T = Tc, when QGP and hadronic matter co-exist, is a 
direct consequence of the first order phase transition. Apart from the role in QGP 
diagnostics, the possibility of the mixed phase affects also the bulk features of the 
evolution process. 

In the mixed phase, the relative proportion of QGP and hadronic matter must 
be a function of time; initially the system consists entirely of QGP and at the end, 
entirely of hot hadronic matter. If we denote the fraction of the QGP by /q(t), 
then the entropy in the mixed phase (smix) can be expressed as, 

Smix = /Q(r)4 + /H(r)s^H (7.21) 

such that at r = tq, fq = 1 and at r = th, fH = ^ — fQ = ^ and the life time 
of the mixed phase Tj™™'^ is th — tq. Here Sq (s^) denotes the entropy density of 
QGP (hadronic) phase at Tc. Since scaling law governing the variation of s(r) must 
continue to hold also in the mixed phase, substituting Eq. ( [7.21|) in Eq. (|7.8|) we 
obtain for Tj > Tc, 

/«m-^(^?-0-,^(t-0 (-^' 

where r (= gq/ges) is the ratio of the degeneracy of QGP phase and the effective 
degeneracy in the hadronic phase. In the above equation we have used the relation 
Th = rrq, obtained as a result of (1 + 1) dimensional isentropic expansion. 

The quantity fq expressed as /q(t) = (s — s'jj)/{sQ — sjj) is the volume fraction 
of the QGP sector in the mixed phase and similarly fn = {sq — s)/{sq — s^) is the 
volume fraction of the hadronic sector in the mixed phase. These quantities (/q(t) 
and fnij)), will be required to evaluate the electromagnetic probes from a evolution 
scenario, QGP— > mixed phase— > hadronic phase— > freeze-out, in the next section. 
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If Ti = Tc, i.e. if the system is formed in the mixed phase with a fraction /o of 
the QGP phase then fg is given by [|^, |142| 



fQir) = — 



1 



;i + (r-l)/o)--l 

r 



(7.23) 



The mixed phase ends at a proper time = {1 + {r — l)/o)rj. In case of Sj < s^, 
the value of fni'^) is always unity. 

To make our discussion more specific, consider Pb + Pb collisions at CERN 
SPS energies. If we assume that the matter is formed in the QGP phase with two 
flavours {u and d), then gk = 37. Taking dNT^/dy = 600 as measured by the NA49 
Collaboration ||143|| for Pb + Pb collisions, we obtain Tj = 185 MeV for = 1 fm/c. 
We have taken Ty = 130 MeV [|144|, p.45|| in our calculations. We also consider central 



collisions of Pb + Pb at the RHIC energies which correspond to about 200 GeV/A 
in the centre of mass system. The particle rapidity density in the central region 
is taken as 1735 for RHIC. The corresponding initial temperature is (by assuming 
that Ti = l and QGP initial state) Ti= 265 MeV . 

8 Results 

8.1 Hadronic properties at non-zero temperature 

In the Walecka model the effective nucleon mass at T 7^ has been evaluated in the 
Relativistic Hartree Approximation (RHA). Then the p and u masses are computed 
by evaluating their self energies due to p — N — N and uj — N — N interactions at 



finite temperature. The following values of the coupling constants and masses ||104 



have been used in our calculations: Kp = 6.1, Qp^f^ = 6.91, mo-= 458 MeV, nip = 770 
MeV, Mn = 939 MeV, c/^^^ = 54.3, = 0, and glj^^ = 102. In Fig. (§) we depict 
the variation of vector meson masses as a function of temperature in the Walecka 
model along with the BR and Nambu scaling scenarios. The parametrized forms 
of the effective masses are given in Eqs. ( 4.5U| ) and ( p.49| ). The variation of mass 



in the Walecka model and BR scaling is slower than the Nambu scaling scenario. 
At higher temperature the Walecka model calculation and the BR scaling (near Tc) 
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Figure 6: Variation of vector meson mass with temperature for BR (long-dashed 
hne), Nambu (dot-dashed hne) scahng with Tc=200 MeV and in the Walecka model 
for p (solid line) and uj (dotted line). 

tend to converge. Such a small difference in the mass variation in the above two 
scenarios may not be visible through the photon spectra. We also note at this point 
that in the Walecka model p and uj masses show different rate of reduction due 
to different values of their coupling constants with the nucleons. 

In Figs. (^ and (|^) the change in the p and uj spectral functions at non-zero 
temperature has been displayed. Here as well as in the following 'free' mass will 
indicate the physical mass in vacuum, e.g. 770 MeV for the p meson. Ap (in units of 
e) is obtained by multiplying Eq. ( |3.19|) by Svr. is obtained analogously. We have 
used the interaction Lagrangians (|3.^), ( |3.10| ) and ( |4.1D for this purpose. The shifts 
in both the spectral functions towards the lower invariant mass region correspond 
to the reduction of their masses due to thermal interactions (see Fig. P). The broad 
UJ peak arises due to its interaction with the thermal pion in the heat bath. The 
reaction vr — >■ vr tt contributes dominantly to the survival probability of the uj in 
the medium. 

Before we proceed further a few comments on Walecka model calculations are in 
order. In this model the major contribution to the medium effects on the p and uj 
mesons arises from the nucleon-loop diagram. For the dressing of internal lines in 
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Figure 7: Spectral function of p meson in the Walecka model. Solid (long dashed) 
line corresponds to T =180 MeV (T =150 MeV). The spectral function in vacuum 
is shown by the dotted line. 




0.3 0.8 1.3 

M (GeV) 



Figure 8: Same as Fig. |^ for the u meson. 
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Figure 9: Spectral function for the isovector channel extracted from e^e~ collisions 
(dotted line) as a function of invariant mass. The dashed (solid) line indicates the 
spectral function when rrip and ujq vary according to BR (Nambu) scaling. 



matter we restrict ourselves to the Mean Field Theory (MFT) to avoid a plethora 
of diagrams and to maintain internal consistency. It has been shown ^ that 
the change in the p mass due to p — tt — vr interaction is negligibly small at non- 
zero temperature and zero baryon density. Therefore the change in the p meson 
mass due to p — vr — vr interaction is neglected here. At finite baryon density, the 
dynamics is more involved due to the medium effects on the p—Ti — n vertex, the pion 
propagator coupled with delta-hole excitation, and the coupling of the p-meson with 
A^*-hole excitations |146|, [M^, [148|, ||, |14|, |15|, |15| (see also the review, Ref.p2|.) 
The major effect of such medium modifications is to broaden the p-peak as well as 
to produce complicated structure around the peak. Since in the present work we 
restrict our calculations within the realm of MFT, i.e the internal nucleon loop in 
the p and u self energy is modified due to tadpole diagram only, the inclusion of 
vertex corrections, modification of the pion propagator and the inclusion of baryon 
resonances are not considered here. Also, in the present work we restrict to zero 
baryon density. 

In Fig. (|) the spectral function (Svr times Eq. ( |6.33| )) for the isovector (p) channel 
is plotted as a function of invariant mass at T = 150 MeV and = 160 MeV. We 
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Figure 10: Same as Fig. | at T = 180 MeV and = 200 MeV. 



find that both the peak and the continuum threshold of the spectral function move 
towards lower invariant mass. However, in case of Nambu scaling the shift is more 
compared to BR scaling. In the Nambu scaling scenario the peak of the spectral 
function and the continuum are not well separated; a merging of the two would take 
place at T = Tc. This could possibly indicate the onset of a deconfinement phase 
transition. Fig. ([10|) shows the spectral function at T = 180 MeV and = 200 
MeV. Due to a larger separation between Tc and T compared to the previous case the 
peaks in the spectral function in all the cases are well separated from the continuum. 

In Figs. (|Tl|) and (|12]) the spectral functions for the isoscalar {uj) channel have 
been depicted. In both the cases the peak in the spectral function is distinctly 
visible in all the mass variation scenarios. The larger width in the isoscalar channel 
is due to the combined processes — Stt and ujtt ^ tt vr as discussed before. 

The spectral functions for the vector mesons both in the isoscalar and isovector 
channels are plotted in Fig. (|T3) at a temperature T ~ Tc. As expected from the 
scaling law the peak has vanished due to its overlap with the continuum. All the 
hadrons in the thermal bath have melted to their fundamental constituents - the 
quarks and gluons. Such a spectral function would indicate a transition from hot 
hadronic matter to QGP. This behaviour should, in principle, be reflected in the 
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Figure 11: Same as Fig. ^ for the isoscalar channel. 
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Figure 12: Same as Fig. ^ for the isoscalar channel. 
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Figure 13: Spectral functions for isovector (p) and isoscalar {u) channels at Tc . 

dilepton spectrum originating from these channels. Such a broad spectral represen- 
tation without any peak may be compared with that of the huge decay width of the 
p meson due to its interaction with the baryonic medium P3| , |153| . 

In Fig. the shift in the pole position of the p-spectral function is depicted 
for the Linear Sigma Model (LSM), Non-Linear Sigma Model (NLSM), and Hidden 
Local Symmetry (HLS) approach. The width of the p meson has been calculated for 
the process p —>■ vr tt at non-zero temperature. For the NLSM and HLS interactions 
the p-mass increases by an amount 90 MeV and 10 MeV respectively. For the 
enhancement of p mass in the NLSM, a larger phase space is available for the decay 
process p ^ -mi consequently the p appears to be broader in this case compared to 
HLS interaction. On the other hand, for the gauged LSM the p mass reduces by 
about 45 MeV at T = 150 MeV. It may be noted, as mentioned in section 5.1, that p 
mass decreases in gauged LSM for low temperatures and increases for temperatures 
in the vicinity of the chiral phase transition. 

8.2 Static photon spectra 

In this section we consider the photon spectra from hot hadronic matter and QGP. 
The medium effects enter through the masses and decay widths of the particles 
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Figure 14: Shift in the pole position of the p spectral function for gauged Linear 
and Non-Linear Sigma Models and Hidden Local Symmetry Lagrangian at T = 150 
MeV. 

participating in the photon producing reactions. It is well known [|lT| (see also 
Refs. 0) that the reactions vrp — > 7r7,7r7r p'j , im rj'y , 111] — > , 
and the decays p ^ nTcy and c<j ^ vr 7 are the most important channels for photon 
production from hadronic matter in the energy regime of our interest. We have also 
included those reactions which produce photon via intermediary axial vector meson 
Oi as discussed earlier. Non-zero width of vector and axial vector mesons in the 
intermediate state has been taken into account. While evaluating the photons from 
QGP we have considered both one loop and two loop contributions to the photon 
self energy as shown in Figs, (g) and (|). 

The total photon emission rate from QGP and hadronic matter at T = 160 
MeV is plotted in Fig. (|T5|) as a function of the energy of the emitted photon for 
different values of strong charge Qs in the QGP phase and for various mass variation 
scenarios in the hadronic sector. The photon production rate from QGP has been 
evaluated in the HTL approximation, which is valid if the hard and soft scales are 
well separated, i.e. for gs « 1 (which corresponds to as « 0.08), QCD lattice 



calculations ||154|| however suggest that as ~ 0.2 — 0.3 at temperatures achievable 
in URHIC. This means that the extrapolation of the results obtained under HTL 
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Figure 15: Thermal photon spectra at T = 160 MeV. Sohd dots (square) indicates 
photon emission rate from QGP with both one loop and two loop contributions as 
evaluated by Kapusta et al and Aurenche et a/ respectively for gg = 2(0.8). Dotdash 
line represents photon spectrum from hot hadronic gas without medium effects. The 
result with the in-medium effects within the scope of the Walecka model calculations 
is shown by long dashed line. Dotted (solid) line indicates photon spectrum with 
BR (Nambu) scaling mass variation scenario. 

approximation to higher values of Qs (or as) corresponding to lattice simulation 
may be dubious. We have evaluated the photon spectra for two values of the strong 
coupling constants Qs = 0.8 (solid square) and 2 (solid dots) to demonstrate the 
sensitivity of the photon spectra to the value of the strong charge and to show the 
uncertainties involved in the problem. In the hadronic sector the photon yield is 
seen to be enhanced compared to the case when the effects of the thermal interaction 
on the hadronic properties are neglected. This is true for almost the entire energy 
range of the emitted photon under consideration. As a result of the similar mass 
shift in the Walecka model and BR scaling the photon spectra in these two scenarios 
(long-dashed and dotted lines respectively) have a negligible difference, whereas the 
enhancement in the spectrum due to hadronic mass shift according to Nambu scaling 
is clearly visible (solid line). 

In Fig. ([16|) we show the photon emission rate at T = 180 MeV. Photon spectra 
from hadronic matter with mass variation according to the Nambu scaling scenario 
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Figure 16: Same as Fig. (pi) at T = 180 MeV and = 2. 



overshine the photons from QGP even for a larger value of Qs (~ 2). 

At this stage one might ask: At a fixed T which one is brighter - the hot hadronic 
gas or QGP?. Well, within the scope of the present work, the answer depends on (i) 
the value of the strong coupling constant, (ii) the degree of hotness of the medium 
and (iii) how adversely the hadrons are affected in the medium. 

In Fig. (|1^) the effect of the form factor on the reaction tttt p'-f has been 
demonstrated. We have taken same monopole form factor for both the vrvrp and 
iriT'j vertices |jTl| to suppress the contribution from very high momentum region 
where the quark structure of the hadrons could be relevant. The Ward-Takahashi 
identity has been used to obtain the dressed propagator. The in-medium mass of the 
p meson has been taken from Walecka model calculations. The form factor effects for 
the above reaction reduces the photon production rate by about 10-15%. In view 
of the experimental uncertainty of the photon spectra measured in URHIC {e.g. 
Ref. ||155|| ) such effects are not relevant at present. Therefore, we have neglected it 
in the following discussions. 

In Fig. (|l^) the emission rate of photon is shown for T = 150 MeV. The change 
in the photon spectra with in-medium masses calculated in the framework of gauged 
LSM, NLSM and HLS Lagrangians are compared with that of vacuum values. An 
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Figure 17: The effect of the monopole form factor on the photon emission rate 
from the reaction tttt pj. 



increase of p mass in the NLSM reduces its number due to Boltzmann suppression, 
which leads to a suppression in the photon emission rate (dotted hne). The pro- 
duction rate is enhanced due to a reduction in the p mass (sohd hne) in LSM. (We 
must remember that the p mass decreases in gauged LSM for low temperature and 
increases for temperatures close to the chiral transition temperature T^. Therefore, 
for T ~ it will show a reduction in photon emission rate and the net yield would 
be a superposition of all temperatures, from initial to freeze-out). The change in the 
mass of p is so small for HLS approach (short-dashed line) that the production rate 
is almost indistinguishable from the spectra with vacuum masses of the hadrons. 
In Fig. (|18D we have also demonstrated how the photon spectra is modified for a 
drastic change in the width of the p meson (Fp ~ 400 MeV, such a large value of 
the width is obtained in Ref. | |153|| although, at non-zero baryon density) without 
any appreciable change in the pole mass {nip ~ 770 MeV). Such drastic broadening 
of the p meson have been proposed in Ref. |Q |153|| . We observe that the effects of 
such modifications in the properties of p on the photon spectra is rather negligible 
(long-dashed line). We will see later that such changes in the properties of p modify 
the dilepton spectra drastically. 
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Figure 18: The change in the photon spectra due to the finite temperature ef- 
fects on the hadronic masses in Linear, Non-Linear Sigma Model and Hidden Local 
Symmetry approach at T= 150 MeV. 

8.3 Static dilepton spectra 

In Fig. ([19|) we display the invariant mass distribution of e~^e~ pair. The dilepton 
yield from qq annihilation is denoted by solid dots. Dotted line indicates the result 
obtained from the parametrization of the electromagnetic current-current correlation 
function in the 'p' and 'cj' channels, when the medium effects are ignored. A large 
shift towards the lower invariant mass region of the p peak is seen in the Nambu 
scaling (solid line) as compared to the BR scaling (dash line) consistent with the 
relative shift in the spectral functions in the two cases as discussed before. In the 
Walecka model calculations the relevant reactions are vrvr — >■ e~^e~, p e^e~ and 
oj e'^e~ (dotdash line) ||132| , |133|| . The two peaks corresponding to p and u masses 



are visible in the spectra. The separation between the two peaks is due to different 
mass shift of the p and u. Measurement of such separation in hadronic masses 
(Am = — m*) would signal the in-medium effects. Validity of such results could 
be tested in URHIC by the CERES [p.56|| collaboration in future. Similar shift at 
zero temperature but finite baryon density could be detected by HADES | |157| | and 
CEB AF . Effects of the continuum on the dilepton spectra is clearly visible for 
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Figure 19: Thermal dilepton spectra at T = 150 MeV. Solid dots indicates dilepton 
emission rate from QGP. Dotted line represents dilepton yield from hot hadronic 
gas without medium effects. The result with the in-medium effects within the scope 
of the Walecka model calculations is shown by the dot-dashed line. Long dashed 
(solid) line indicates dilepton spectrum with BR (Nambu) scaling mass variation 
scenario. 



M > 1 GeV (please note that the value of the continuum threshold in vacuum is 1.3 
GeV). Due to the continuum contribution the dilepton rates from hadronic matter 
and QGP shine equally brightly in the mass range M > 1 GeV. The lepton pair 
spectra at T = 180 MeV is shown in Fig. (PUj). Since the effective mass of the p in 
the Walecka and BR scaling scenario is almost same in this case (see Fig. @), the 
corresponding rates are very similar near the p peak. 

The dilepton invariant mass distribution at T = Tc is shown in Fig. (0). All the 
peaks in the spectrum have disappeared as expected. The rates obtained from the 
electromagnetic current-current correlator is close to the rate from qq annihilation, 
indicating that the qq interaction in the vector channel has become very weak, 
signaling the onset of deconfinement |[L58| , |159 1 . 

In Fig. (^) we compare the dilepton emission rate at T = 150 MeV for vacuum 
mass of p with the rates where the effective masses are obtained in the framework 
of gauged LSM, NLSM and HLS approach. The positive shift of p mass in NLSM 
is reflected in the peak position of the spectra towards larger value of M (dotted 
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Figure 20: Same as Fig. (|T7D at T = 180 MeV. 




Figure 21: Same as Fig. at T = T, MeV. 



81 




0.5 1 1.5 

M (GeV) 



Figure 22: The change in the invariant mass distribution of lepton pairs due to 
the finite temperature effects on the hadronic masses in Gauged Linear, Non-Linear 
Sigma Model and Hidden Local Symmetry approach at T= 150 MeV. The long- 
dashed hne indicates dilepton spectra for Fp = 400 MeV and rup — 770 MeV (see 
text). 

line). A very small change in the p mass in HLS approach does not cause any visible 
change in the invariant mass distribution of dilepton. The long-dashed line indicates 
lepton pair distribution for rUp = 770 MeV and Tp = 400 MeV. The large width of p 
leads to the disappearance of the p peak from the spectra, indicating that p ceases 
to exist as a quasi-particle. However, as mentioned before, it is interesting to note 
that the photon spectra is "insensitive" to such drastic broadening of the p meson. 

8.4 Photon and dilepton spectra with space-time evolution 

As mentioned before the basic aim of URHIC is to distinguish between the two 
possibilities: 

A -|-A ^QGP^Mixed Phase^Hadronic Phase 

or 

A -|- A ^Hadronic Phase 

The former (latter) case where the initial state is formed in QGP (hadronic) 
phase will be called the 'QGP scenario' ('no phase transition scenario'). In the 
following we will compare the photon and dilepton spectra originating from these 
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two scenarios with and without medium effects. 

The observed photon and dilepton spectra originating from an expanding QGP or 
hadronic matter is obtained by convoluting the static (fixed temperature) rate with 
the expansion dynamics. The basic ingredients required for a system undergoing 
rapid expansion from its initial formation stage to the final freeze-out stage with 
or without phase transition have been discussed in section 6. For the QGP sector 
we use a simple bag model equation of state (EOS) with two flavour degrees of 
freedom. The temperature in the QGP phase evolves according to Bjorken scaling 
law T = TfTi. The cooling law in the hadronic sector is quite different from 
that of the QGP because of the presence of massive hadrons. These hadrons redress 
themselves in the medium thereby changing their vacuum masses. This phenomenon 
must be taken into account in the evolution dynamics through the equation of state. 
We do this by introducing temperature dependence in the statistical degeneracy 
which takes care of the mass varying with temperature. 

In Fig. (^) we depict the variation of effective degeneracy as a function of 
temperature with and without medium effects on the hadronic masses for various 
scenarios. We observe that for T > 140 MeV the effective degeneracy becomes 
larger due to the reduction in temperature dependent masses compared to the free 
hadronic masses. Physically this means that the number of hadrons in a thermal 
bath at a temperature T is more when in-medium mass reduction is taken into 
account. Eq. (|7.2CI|) implies that for a given pion multiplicity the initial temperature 
of the system will be lower (higher) when medium effects on hadronic masses are 
considered (ignored). This is clearly demonstrated in Fig. (|2^) where we show the 
variation of temperature with proper time for different initial conditions. The solid 
dots indicate the scenario where QGP is formed initially at Tj = 185 MeV and cools 
down according to Bjorken law upto a temperature Tc at proper time tq, at which 
a phase transition takes place; it remains constant at Tc up to a time th = 9.4 fm/c 
after which the temperature decreases as T = 0.247/r°'^^'^ (when medium effects 
are taken from Walecka model) to a temperature Tf{= 130 MeV). If the system 
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Figure 23: Variation of effective degeneracy as a function of temperature. 

is considered to be formed in the hadronic phase then the initial temperature is 
obtained as Tj = 220 MeV (270 MeV) when in-medium effects on the hadronic 
masses from Walecka model is taken into account (ignored). The corresponding 
cooling laws are displayed in Fig. (p^. The above parametrizations of the cooling 



law in the hadronic phase have been obtained by solving Eq. ( |7.11| ) self consistently. 
An initial state with the vanishing meson masses at Tj = 195 MeV (xj = 1 fm/c) 
could be realised in the case of BR and Nambu scaling scenarios for the value of 
pion multiplicity, dN/dy = 600. 

In Table 1 we quote the values of the initial temperatures obtained by assuming 
various mass variation scenarios. The value of initial thermalization time has been 
assumed as 1 fm/c both for SPS {dN/dy = 600) and RHIC {dN/dy = 1735) energies. 
tq {th) indicates the starting (end) point of the mixed phase, th — tq is the life 
time of the mixed phase in a first order phase transition scenario. and 6 dictate 
the variation of temperature with proper time for the hadronic matter according to 
the cooling law T = ^O/t^ . The values of 5 indicate a slower cooling in the hadronic 
phase as compared to that in the QGP phase (T ~ l/r° '^^). 
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Figure 24: Variation of temperature as a function of proper time. 
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Ti (MeV) 


Th (fm) 




6 


th (fm) 




6 


free mass 


270 


10.8 


0.267 


0.215 


31.9 


0.337 


0.215 


Walecka 


220 


9.4 


0.247 


0.194 


27.6 


0.305 


0.194 


BR 


195 


8.2 


0.236 


0.184 


23.9 


0.288 


0.185 


Nambu 


195 


4.7 


0.203 


0.151 


13.9 


0.239 


0.152 



Table 1 : Values of initial temperatures and various time scales for SPS and RHIC 

energies. 

Having obtained the finite temperature effects on hadronic properties and the cooling 
laws we now integrate the rates obtained in the previous sections over the space-time 
evolution of the collision. We must account for the fact that the thermal rates are 
evaluated in the rest frame of the emitting matter and hence the momenta of the 
emitted photons or dileptons are expressed in that frame. Accordingly, the integral 
over the expanding matter is of the form 

freeze — out 



dN 



dV 



formation 

where dV stands for invariant phase space elements: d^p/E for photons and d^q 
for dileptons. E* is the energy of the photon or lepton pair in the rest frame of 
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the emitting matter and T{x) is the local temperature. In a fixed frame like the 
laboratory or the centre of mass frame, where the 4- momentum of the photon or 
lepton pair is = {E, q) and the emitting matter element d^x moves with a velocity 

= 7(1, v), the energy in the rest frame of the fiuid element is given by E* = u^q^. 

In a first order phase transition scenario the photon and dilepton spectra from 
a (1 + 1) dimensionally expanding system is obtained as 



dN 



'dB\ 



QGP 



+ 



[dr) 

dR\ s 



S — S 



H 



QGP ^ 



H 



6(4 - .)e(. - s'^h) 

T dr drj 



12) 



where Ra is the radius of the nuclei and 9 functions are introduced to get the 
contribution from individual phases, and Sq are defined in section 7.2. 

As discussed earlier, g^s is obtained as a function of T by solving Eq. (|7.19|) . A 
smaller (larger) value of Qes is obtained in the free (effective) mass scenario. As a 
result we get a larger (smaller) initial temperature by solving Eq. ( |7.2(]| ) in the free 
(dropping) mass scenario for a given multiplicity. Naively we expect that at a given 
temperature if a meson mass drops its Boltzmann factor will be enhanced and more 



of those mesons will be produced leading to more photons p3| , |160|| . However, a 
larger drop in the hadronic masses results in smaller initial temperature, implying 
that the space time integrated spectra crucially depends on these two competitive 
factors. Therefore, with (without) medium effects one integrates an enhanced (de- 
pleted) static rate over smaller (larger) temperature range for a fixed freeze-out 
temperature (Tj = 130 MeV in the present case). In the present calculation the 
enhancement in the photon emission due to the higher initial temperature in the 
free mass scenario (where static rate is smaller) overwhelms the enhancement of the 
rate due to negative shift in the vector meson masses (where the initial temperature 
is smaller). Accordingly, in the case of free mass ( Nambu scaling) scenario the 
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Figure 25: Total thermal photon yield corresponding to dN/dy = 600 and Tj = 1 
fm/c. The solid (long-dash) line indicates photon spectra when hadronic matter 
formed in the initial state at Tj = 195 MeV (Tj = 220 MeV) and the medium effects 
are taken from Nambu scaling (Walecka model). The dotted line represents the 
photon spectra without medium effects with Tj = 270 MeV. The solid (dotted) line 
with solid dots represent the yield for the 'QGP scenario' when the hadronic mass 
variations are taken from Nambu scaling (free mass). 

photon yield is the highest (lowest). In case of the Walecka model, the photon yield 
lies between the above two limits. This is demonstrated in Fig. (p5|). 

In the 'QGP scenario' the photon yield with in-medium mass is lower than 
the case where free masses of hadrons are considered. However, the difference is 
considerably less than the 'no phase transition scenario'. This is because, in this 
case the initial temperature is determined by the quark and gluon degrees of freedom 
and the main difference between the two is due to the different lifetimes of the mixed 
phase. In Fig. (pSf), the photon spectra from 'QGP scenario' is compared with that 
from 'no phase transition scenario'; the latter overshines the former. 

The space time integrated dilepton spectra for the 'QGP scenario' and 'no phase 
transition scenario' with different mass variation are shown in Fig.(p6D. The shifts 
in the invariant mass distribution of the spectra due to the reduction in the hadronic 
masses according to different models are distinctly visible. Similar to the photon 
spectra, the dilepton spectra from 'no phase transition scenario' dominates over the 
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Figure 26: Total thermal dilepton yield corresponding to dN/dy = 600 and 7^ = 1 
fm/c. The solid (long-dash) line indicates dilepton spectra when hadronic matter 
formed in the initial state at Tj = 195 MeV (Tj = 220 MeV) and the medium 
effects are taken from Nambu scaling (Walecka model). The dotted line represents 
the spectra without medium effects with Tj = 270 MeV. The square (solid dots) 
represent the yield for the 'QGP scenario' when the hadronic mass variations are 
taken from Nambu scaling (Walecka model). 

'QGP' scenario for invariant mass beyond p peak. 

The space time integrated photon yield corresponding to the static (fixed tem- 
perature) rate shown in Fig. (|l^) (when in-medium effects are taken from gauged 
Linear and Non-Linear Sigma models and Hidden Local Symmetry) are not dis- 
played in Fig.(p5D because the resulting spectra does not differ appreciably from 
those already shown and it makes the Fig. (p5l) clumsy. For the same reason the 
corresponding dilepton spectra are also not shown. 

Finally we study the electromagnetic probes for RHIC energies. At RHIC a 
scenario of a pure hot hadronic system within the format of the model used here, 
appears to be unrealistic. The initial temperature considering free hadronic masses 
turns out to be ~ 340 MeV whereas for the other extreme case of massless hadrons 
it is ~ 290 MeV. With temperature dependent masses the initial temperature will 
lie somewhere between these two values. For such high temperatures, clearly a hot 
dense hadronic system cannot be a reality, the hadrons would have melted away 
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Figure 27: Variation of temperature as a function of proper time. The initial tem- 
perature has been determined by assuming 'QGP scenario'. The initial temperature 
Ti = 265 MeV for n = 1 fm/c and dN/dy = 1735. 

even for lower temperatures. Thus, for RHIC we have treated the case of a QGP 
initial state only. The temperature profile for RHIC is depicted in Fig. (^) where 
we observe that the length of the plateau, which indicates the life time of the mixed 
phase r^{^ = th — tq, depends on the masses of the hadrons in the hadronic phase. 
The effective degeneracy plays an important role here. At the transition point there 
is a large decrease in the entropy density. This decrease has to be compensated by 
the expansion (increasing the volume) to keep the total entropy constant. Since we 
are considering (1+1) dimensional expansion this change in the entropy density will 
be compensated by increasing r (sr = const.). We have seen earlier (Fig. (PBD) that 
the effective degeneracy in the hadronic phase is the largest for the Nambu scaling 
and smallest for the free mass scenario, resulting in smallest (largest) discontinuity 
in the entropy density for the former (latter) case. Consequently the time taken for 
the system to compensate the decrease of the entropy density in the Nambu scaling 
scenario is smaller as compared to free mass case. Hence the life time of the mixed 
phase for the Nambu scaling case is smaller than all other cases. 

The thermal photon spectra for RHIC is displayed in Fig. (pSj) . The solid line 
represents the total thermal photon yield originating from initial QGP state, mixed 
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Figure 28: Thermal photon spectra at RHIC energies. 

phase and the pure hadronic phase. The short dash hne indicates photons from 
quark matter (QM) (= pure QGP phase + QGP part of the mixed phase) and the 
long dash line represents photons from hadronic matter (HM) (= hadronic part of 
the mixed phase + pure hadronic phase). In all these cases the effective masses of 
the hadrons have been taken from Nambu scaling. For pt > 2 GeV photons from 
QM overshines those from HM since most of these high photons originate from 
the high temperature QGP phase. We arrive at the similar conclusions when in- 
medium effects in the Walecka model is considered (not shown in Fig. (|28|) to avoid 
clumsiness). The dotted and the dotdash lines indicate photon yields from QM and 
HM respectively with free masses in the hadronic sector. The HM contribution for 
the free mass is larger than the effective mass (Nambu) scenario because of the larger 
value of the life time of the mixed phase in the earlier case (see Table 1). We note 
here that for p^ > 2 GeV, photons from QM overshine those from HM irrespective 
of the models used for calculating the in-medium modifications of the hadrons. 

Thermal dilepton yield at RHIC energies for QGP initial state and for different 
mass variation scenarios are shown in Fig. (pQ]). The shape of the peak in the dilepton 
spectra for Walecka model is slightly different (broader) from that of Nambu scenario 
because of the larger mass separation between p and u mesons in the former case 
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Figure 29: Thermal dilepton spectra at RHIC energies. 

(see Fig. (^). The shift in the peak of the dilepton spectra towards lower invariant 
mass for both the Walecka model and Nambu scaling scenario is clearly visible. 
However, the dilepton yield from QGP dominates over the contribution from the 
hadronic phase for M > 1.4 (M > 2.3) GeV for Walecka model (Nambu scaling). 

9 Summary and Outlook 

In the present work we have reviewed the formulation of the production of photon 
and lepton pair from QGP and hot hadronic gas based on finite temperature field 
theory. The changes in the spectral functions of the hadrons appearing in the 
internal loop of the photon self energy diagram have been considered in the Walecka 
model, gauged linear sigma model, non-linear sigma model, hidden local symmetry 
approach and QCD sum rule approach. The hadronic spectral functions (in vacuum) 
for the isovector and isoscalar channel have been constrained from the experimental 
data of e^e" hadrons. Due to the lack of our understanding of the critical 
behaviour of scalar and tensor condensates we have parametrized the vector meson 
masses and continuum threshold as a function of temperature according to BR and 
Nambu scaling. 

We observe that the in-medium effects on the hadronic properties within the 
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frame work of the Gauged Linear and Non-Linear Sigma Model, Hidden Local Sym- 
metry approach are too small to affect the electromagnetic spectra substantially. 
However, the shift in the hadronic properties of different magnitude within the 
frame work of the Walecka model, Brown-Rho scaling and Nambu scaling scenarios 
are prominently visible through the low invariant mass distribution of dilepton and 
transverse momentum spectra of photon at SPS energies. The photons from 'no 
phase transition scenario' outshine those originating from the 'QGP scenario' for 
the entire range of px at SPS. At RHIC energies a scenario of a pure hot hadronic 
system appears to be unrealistic because of the very high initial temperature ob- 
tained within the format of the model used in the present work. We observe that at 
RHIC energies the thermal photon (dilepton) spectra originating from Quark Gluon 
Plasma overshines those from hadronic matter for high transverse momentum (in- 
variant mass) irrespective of the models used for evaluating the finite temperature 
effects on the hadronic properties. 

We note that Walecka model calculation gives different mass shift for p and uo 
mesons (because p and uo couples to nucleons with different strength). The disen- 
tanglement of the p and uo peaks in the dilepton spectrum resulting from URHIC 



would be an excellent evidence of in-medium mass shift of vector mesons 156 



and/or validity of such model calculations for the situation under consideration. It 
is also interesting to note that the dilepton spectra is affected both by the changes 
in the decay width as well as in the mass of the vector mesons. However, the photon 
spectra is affected only by the change in the mass of the vector mesons but rather 
insensitive to the change in its width. The effects of the continuum on the dilepton 
spectra are seen to be substantial. 

In the following we would like to reiterate some of the assumptions made in the 
present calculations. Further work is needed for the justifications of these assump- 
tions. 

The exact value of the critical temperature (T^) for the deconfinement phase 
transition is still uncertain. However, recent lattice simulation [§] for two flavour 
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QCD indicates a value of Tc for chiral transition ~ 130 — 160 MeV. We have taken 
Tc = 160 MeV, although till now it is not known whether the values of for the 
chiral and deconfinement transition are the same or not. The value of the initial 
thermalization time is unfortunately also an unknown quantity (though consid- 
erable progress has been made for the understanding of the initial condition ||161|] , 
a lot is yet to be done). We take = 1 fm/c as a canonical value following 
Bjorken ||135|| . A similar value of Tj has been considered in the literature, e.g. see 
Refs. m, |14|, |T6|, |16| . 

The photon production from QGP has been evaluated using HTL resummation 
based on the assumption « 1, which is unlikely to meet in URHIC even at the 
highest energy to be available at the CERN LHC in future. The strong coupling 
constant is likely to attain a value ~ 2 at RHIC/LHC. Evaluation of the photon 
spectra at such high values is a formidable task. In this respect the development of 
methods suitable for addressing non-perturbative effects near and above the QCD 
phase transition point is of paramount importance. Extension of the self-consistent 
resummation scheme developed in 0^ theory ||163|| to non-abelian gauge theory ||164 



165|| would be a very important step towards the understanding of the phenomena 



near the QCD phase transition. 

In this work we have neglected transverse expansion of the evolving matter. This 
is because our main emphasis has been to study the shift in the photon and dilepton 
spectra due to the change in the hadronic spectral function at non-zero temperature. 
Inclusion of the transverse expansion will shift the momentum distribution of the 
photon both for the free mass and in-medium mass scenarios but the relative shift 
will remain approximately unchanged and it is this relative change in which we are 
interested. The effects of the non-zero baryonic chemical potential are ignored in 
the present work, calculation addressing these issues is in progress [|166| . 



Throughout this work we have assumed thermal equilibrium, which may not 
be realized practically | |167| , |168| , |169| , |170| , |171| |. Unfortunately, although consid- 



erable progress has been made [p.72| , |173| , |174|| , the general techniques for solving 
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non-equilibrium quantum field theoretical problems is still in the early stages of 
development |[175| . 

In spite of some of the remaining uncertainties, the progress in this field is 
remarkable, especially when one considers that not a tremendous number of experi- 
ments to seek out QGP have been performed. Many of the ambiguities pertaining to 
the "pre-data" theory have been removed by the experimental data from AGS and 
SPS. We have all the reasons to look forward to RHIC and LHC where the initial 
energy density will be so high that the formation of QGP is almost inevitable. 
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Appendix - Thermal Propagators 

In this appendix we will briefly discuss the in-medium (thermal) propagators in 
Quantum Field Theory p.76| , |177|| which have been used extensively in this 



work. We will begin by first defining the propagators in vacuum. 

The free propagator of a complex scalar field propagating with a momentum 
p in vacuum is defined as (see e.g. ||178|| ) 

ih\p) = j d*xe'-(T{4,(x)4,'(mo 

(1.1) 



— w? + ie 

The operator T appearing within angular brackets ensures that the field operators 
are time-ordered and the subscript '0' indicates that there are no interactions. The 
vacuum propagator for fermions is defined as 

j4. 



(1.2) 



^G%{p) = J d^xe'P-^T{M^)MO)})o 



p2 _ ^2 _j_ 



where a and (3 denote the spinor indices of the fermion field ip. For massive vector 
particles we define the free propagator in vacuum as 

tDlip) ^ f d'xe'P-^T{A,{x)A,m)o 



tj-g^-'+p^^p-'/m') 

pZ _ _|_ 

where is the massive vector field. For charged vector particles the fields A^ will 
also have isospin indices. It must be noted that the fields 0, ip and A^ appearing 
above are free fields and the expectation values are calculated between noninteract- 
ing vacuum states. The propagators defined through Eqs. (pTf), ( p2D and (|L3|) are 
referred to as Feynman propagators. 

In the presence of interactions these propagators have to be redefined with in- 
teracting Heisenberg fields in place of the free fields and interacting vacua instead 
of the free vacua. The interacting propagator can be expressed in terms of the 
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free(non-mteracting) propagator through perturbation theory. In the scalar case, 
for example, the exact propagator in the presence of interactions. A, is obtained as 



A = Ao + AqIIA 



(1.4) 



Where, 11 is the self energy of the particle due to interactions. This equation is 
known as the Dyson-Schwinger equation for propagators. 

Let us now study the situation in a medium at finite temperature (and density). 
We will be interested in a system in thermal equilibrium. Hence we will assume that 
the interaction slowly switches off as we go into the remote past and the fields become 
noninteracting fields satisfying the free equations of motion. These fields appear in 
the definition of the free propagators in the medium. The thermal propagator has 
more structure than the vacuum one arising due to different combinations of time- 
In the real time formalism there 



ordering on the real time contour |^9|, Q |6( 
are four non-trivial propagator structures possible which are collected in a 2x2 
matrix ||179| , |18CI|| . For scalars the free thermal propagator is defined as 



J d^xe'P-''{T{(j){x)(j)*{0)})^T I ci^a;e^f^({0*(O)0(x)})O 
/ rfSe*f^({0(x)0*(O)})§. / rf^a;e^f^(T{0(x)0*(O)})?. 



(1.5) 



where the operator T denotes anti-time-ordered product. The subscript 'T' indicates 
that a thermal average is being performed. 

In order to obtain the thermal propagators in momentum space one follows 
the usual procedure of expanding the field operators in terms of the creation and 
annihilation operators and making use of the commutation relations between them. 
The four components are then obtained as 



2m5{p — m )rj{p.u) 

<^ — + It 

-2T\%b{j? — m'^)[r](p.u) + 9{—p.u)] 



-2T{i5ijP' — m'^)[rj{p.u) + 6{p.u)] 

: 2'Ki5{p^ — m^)r](p.u) 



(1.6) 



96 



where T]{p.u) = 6{p.u)fBE{z) + 6{-p.u)fBE{-z). Jbe = [e"" - l]'^ is the Bose 
distribution with z = {p ■ u — fi)/T and is the four velocity of the thermal bath. 
We observe that the elements of the matrix propagator Aq are not independent. 
From their definitions one can see that Aj^ and Aq^ can be expressed in terms of Aq 
and Aq^. Also, the Kubo-Martin-Schwinger periodicity condition yields Aq^(p) = 



12 



-Al\p). 

Similarly, for fermions the corresponding thermal propagators are. 



Explicitly, the four components are 



^O(ll) 



ip) 



-J d''xe'P-^{{MO)Mx)})T 
J d^xe^P-^{T{Mx)M^)})T 



+ 2'Ki5{p^ — m^)ri{p.u) 



(1.7) 



gT\p) 

^0(21) 



2Txi{j) + vrL)af3^{j? — m^) [r](p.u) — 9(—p.u)] 
2Tii{jl> + m)ai35{p^ — m^) [rj{p.u) — 6{p.u)] 



ii) + ni] 



0/3 



p2 _ ^2 



It 



+ 2niS{p'^ — m'^)r](p.u) 



{U 



where, rj{p.u) = 9{p.u)fFD{z) + 6{—p.u)fpD{—z), fpo = [e^ + 1] ^, the Fermi- 
Dirac distribution. For the fermions the KMS anti-periodicity condition leads to 



,^12 



21 



Lastly, we define the finite temperature propagators for massive vector particles: 



/ d^xe'P--{T{A^{x)AM])T I d^xe^P--{{A,{0)A^{x)}fT 
J d^xe^P--{{A^{x)A,{0)})^T I d^xe^P-^{f{A^{x)A,iO)})^T 



(1.9) 



where the isospin indices relevant for charged vector particles like the p meson have 
been suppressed. 

It is important to note that the real time propagators as given by Eqs. (I^) and 
(P^) consist of two parts - one corresponding to the vacuum, describing the exchange 
of virtual particle and the other, the temperature dependent part, describing the 
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participation of real (on-shell) particles present in the thermal bath in the emission 
and absorption processes. The temperature dependent part does not change the 
ultra-violet behaviour of the theory as it contains on-shell contributions which has 
a natural cut-off due to the Boltzmann factor. Therefore, the zero temperature 
counter term is adequate for the renormalization of the theory. However, the infra- 
red problem becomes more severe at finite temperature [|^, |6^ . 

Thermal field theory in the real time approach can be reformulated by diago- 
nalising the 2x2 matrix propagators described above. A well-known possibility is 
to diagonalise to a matrix constructed from the Feynman propagators. The non- 
interacting propagator defined by Eq. ( |1.5D can be written as 



Ao = U 



Ao 






-A* 

^0 



U 



:i.io) 



where 



U 



Tj + 9{p ■ u) 



T] + 9{—p ■ u) 



The exact propagators in the medium can be defined analogously as Eqs. (p5|) , (p^ ) 



and ( |I.9|) with interacting Heisenberg fields instead of the free fields. In this case we 
write 



U 



A 






-A* 



U 



where A is the matrix of interacting thermal propagators. Using thermal perturba- 
tion theory A can be expanded in terms of Aq. One obtains 



A = Ao + AqIIA 



(1.12) 



where 11 now is a matrix of self energies 

n = 



n 
-n* 



The diagonal component of Eq. ( [1.1 2| ) is given by 



A = Ao + AqITA 



:i.l3) 



(1.14) 
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In the following we will discuss the vector (spin 1) propagator in some detail. It 
is very similar to the scalar case except now we have to take into account the Lorentz 
structure of the propagator and the self energy. The exact propagator (matrix) D 
can be diagonalised as above and the diagonal element satisfies Dyson equation 

D,, = D% + Dl^W'^D„, (1.15) 

which gives 

D-l = {D%)-'-Yi,,, (1.16) 

where D^^^ is the vacuum Feynman propagator for massive vector particles as defined 
above and Il^y is the self energy which can be written as a sum of two contributions: 

^"'^ = ^Zc + ^Ze,. (1-17) 

where 



nc = (^7"'^-^)nvac(/), (1.18) 

is the vacuum contribution to the self energy, = {uj,p) is the four-momentum of 
the propagating particle. In a thermal bath moving with four-velocity u^, 11^^^^ has 



transverse and longitudinal components EM: 



^Zdi^^P) = A'"'UT,med + B^''UL,raed. (1.19) 

A'^'^ and B^'^ are the transverse and longitudinal projection tensors given by 
1 



(/ _ cu2)(^'^" - u^u"") - p^p'' - ljW + uiu^'p'' + p^'u") 



and 

1 



B 



p2 (jjl _ ^Jl 

which satisfy the following algebra: 



(1.20) 
(1.21) 



B,,B^' = b; 

A.pB'"' = 

p2 
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Using Eqs. ( [i.l6[ - [i.2^ ) the effective propagator becomes 

p) ^^"^ -^M^ |_ PfJ-Pf /T r)o\ 

J^flU — r, OiTT 9 9iTT~'~4.' K^-'^'^} 

where 

nr(L) = ^T(L),raed + (1-24) 

and m, we recall, is the bare mass of the particle. The real part of TIt^l) affects 
the dispersion relation of the particle in the medium. The displaced pole position 
of the effective propagator in the rest frame of the propagating particle (i.e. where 
the three momentum of the particle is zero) gives the effective mass of the particle 
in the medium. The imaginary part of T{t{l) is connected to the decay width. 

A different scheme in the formulation of finite temperature field theory known 



as the 'R/A' formalism ||181|| (see also |p.82|| ), is to diagonalise to a matrix composed 



of retarded and advanced propagators which are known to have better analyticity 
properties than the Feynman ones. In this case the analogue of Eq. (|I.10|) is 



Ao = V 



A? 
A(f 



W. (1.25) 



The matrices V and W depend on the momentum as well as the thermal factor 
containing the distribution functions. Their exact forms are given in Ref. ||181|j . 
The retarded and advanced propagators are defined as 

2A« = j rf^xe^^-^^(xo)([0(x),0(O)])o 

zA^ = / rf^xe*P-^^(-Xo)([0(x),0(O)])o. (1.26) 



For the case of vectors one arrives at the following expression for the effective re- 
tarded propagator at finite temperature: 

where 11;^ and 11^ are respectively the retarded transverse and longitudinal compo- 
nents of the self energy. 

Before we end our discussion on finite temperature propagators let us briefly 
mention about the imaginary time formalism or Matsubara formalism which has 
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been used extensively in the literature. In the imaginary time formalism, the form 
of the propagator at finite temperature remains same as that of vacuum but the 
time component of the four-momentum takes discrete values, i.e. po = 2mTiT{= 
(2n + l)7rzT) for bosons (fermions) with n = — oo to +00, the vertices are the 
same as the zero temperature theory and the loop integral / (i^p/(27r)^ is replaced 
by the sum iTJ2nI d^p/i^Tr)^. There are standard tricks to evaluate the sum over 
the frequencies [^. Another method (known as S ACL AY method) which uses the 
mixed representation of the propagator i.e. it depends on the three momentum 



and Euclidean time has also been used extensively in the literature [Tj], a detailed 
discussion of which can be found in Ref. |183| . The propagators in the imaginary 
time formalism can also be obtained by proper analytic continuation of the real time 



propagators ||184| , |185|| . 
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